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This  dissertation  Is  concerned  with  the  dynamical  analysis  of  an 
elastic  bar  whose  stress-strain  relation  Is  not  monotonic.  Sufficiently 
large  applied  loads  then  require  the  strain  to  jump  from  one  ascending 
branch  of  the  stress-strain  curve  to  another  such  branch.  For  a  special 
class  of  these  materials,  a  nonlinear  Initial -boundary  value  problem  In 
one-dimensional  elasticity  Is  considered  for  a  semi-infinite  bar  whose 
end  Is  subjected  to  either  a  monotonlcally  Increasing  prescribed  traction 
or  a  monotonlcally  Increasing  prescribed  displacement.  If  the  stress  at 
the  end  of  the  bar  exceeds  the  value  of  the  stress  at  any  local  maximum  of 
the  stress-strain  curve  a  strain  discontinuity  or  "phase  boundary"  emerges 
at  the  end  of  the  bar  and  subsequently  propagates  Into  the  Interior.  For 
classically  smooth  solutions  away  from  the  phase  boundary,  the  problem  Is 
reducible  to  a  pair  of  differential -del ay  equations  for  the  two  unknown  func 
tlons  of  a  single  variable.  The  first  of  these  two  functions  gives  the 
location  of  the  phase  boundary,  while  the  second  characterizes  the  dynami¬ 
cal  fields  In  the  high-strain  phase  of  the  material.  In  these  equations 
the  fbrmer  function  occurs  In  the  argument  of  the  latter>v$o  that  the  delays 
In  the  functional  equations  are  unknown..  A  short-time  analysis  of  this 

system  provides  an  asymptotic  description  of  the  emergence  and  Initial  \ 
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♦The  results  communicated  In  this  paper  were  obtained  In  the  course  of  an 
Investigation  supported  In  part  by  Contract  N00014-75-C-0196  with  the  Office 
of  Naval  Research  In  Washington,  D.C. 


propagation  of  the  phase  boundary.  For  large- times,  a  different  analysis 
Indicates  that  the  phase  boundary  velocity  approaches  a  constant  which 
j  depends  on  material  properties  and  on  the  ultimate  level  reached  by  the 
applied  load  as  weTV^  Higher  order  corrections  depend  on  the  detailed 
way  In  which  the  load  Is  appVtod* >  Estimates  for  the  various  dynamical 
field  quantities  are  given  and  a  priori  conditions  are  deduced  which 
determine  whether  the  phase  boundary  eventually  becomes  the  leading 
disturbance. 
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INTRODUCTION 


The  phenomenon  of  phase  changes  has  lately  received  Increased  atten¬ 
tion  In  continuum  mechanics.  In  the  setting  of  one-dimensional  finite 
elasticity,  for  example,  the  modelling  of  phase  changes  Involves  a  strain 
energy  density  which  Is  a  non-convex  function  of  the  displacement  gradient. 
This  leads  to  a  stress-strain  relation  which  Is  not  monotonic.  Consequent 
ly  the  values  of  strain  where  the  stress  response  changes  monotonicity 
partition  the  strain  axis  Into  distinct  Intervals.  We  associate  each  such 
strain  Interval  with  a  distinct  material  phase.  Under  certain  loadings  of 
a  body  composed  of  such  a  material,  the  possibility  exists  that  surfaces 
emerge  across  which  the  strain  Is  dlscontlnous  and  which  separate  regions  of 
different  phase.  The  present  study  Is  motivated  by  questions  raised  by 
the  recent  literature  Involving  solids  with  such  multiple  phases.  These 
Issues  are  addressed  In  detail  In  the  following  section. 

In  subsequent  sections  we  formulate  an  Initial -boundary-value  problem 
associated  with  the  dynamics  of  a  homogeneous,  elastic  bar  whose  strain 
energy  density  Is  such  a  non-convex  function  of  the  displacement  gradient; 
the  objective  Is  to  Illuminate  the  manner  In  which  the  phase  boundaries 
are  generated  and  propagated  In  response  to  given  loading  conditions. 

In  order  to  avoid  the  consideration  of  reflected  waves,  we  consider  a 
semi-infinite  bar  which  we  take  to  be  Initially  undeformed  and  at  rest. 

4 

Two  problems  are  treated:  In  one,  the  traction  at  the  end  of  the  bar  Is 
prescribed  as  a  function  of  time,  corresponding  to  a  load-controlled 
(or  "soft")  dynamical  testing  device;  In  the  second -or  "hard" -dynamical 
problem,  the  displacement  history  of  the  end  of  the  bar  Is  specified. 
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Condltlons  are  then  put  on  the  strain  energy  density  and  on  the  loadings 
so  as  to  focus  attention  on  the  Issue  of  phase  boundary  propagation  with 
a  minimum  of  technical  complication.  We  show  how  certain  loadings  require 
the  emergence  of  a  phase  boundary  at  the  end  of  the  bar  which  subsequently 
propagates  Into  the  Interior. 

For  a  subclass  of  the  materials  considered,  the  location  of  the 
phase  boundary  Is  governed  by  a  pair  of  differential -delay  equations, 
where  the  delay  Is  unknown.  We  first  present  a  short-time  analysis  of 
this  system  In  order  to  describe  the  emergence  and  Initial  propagation 
of  the  phase  boundary.  We  then  study  this  system  at  large  times  and 
show  that— for  the  problems  considered —the  phase  boundary  propagates  at 
a  velocity  which  lies  between  the  sound-speeds  of  the  "low-strain"  and 
"high-strain"  phases  of  the  material.  Here  the  phase  boundary  may  or 
may  not  eventually  become  the  leading  disturbance  to  travel  down  the  bar. 
In  either  case  the  speed  of  propagation  of  the  phase  boundary  approaches 
a  constant  which  depends  on  the  ultimate  level  reached  by  the  applied 
load,  as  well  as  on  the  material. 

All  of  the  analysis  described  above  Is  carried  out  in  detail  for  the 
"soft"  device.  A  summary  of  results  for  both  the  "hard"  and  "soft"  device 
is  presented  In  the  concluding  section.  There  we  also  describe  the 
asymptotic  nature  of  the  associated  dynamical  fields  for  both  short  and 
large  times,  as  well  as  discuss  other  pertinent  aspects  of  the  problem. 


1.  Background 

For  the  one-dimensional  homogeneous  elastic  bar  In  which  the  strain 
energy  density  Is  not  a  convex  function  of  the  strain,  the  stress  will 
not  be  a  monotone  function  of  strain.  Erlcksen  [1]  has  investigated 
equilibrium  configurations  of  bars  whose  stress-strain  relation  is  of 
the  type  depicted  In  Fig.  1,  In  which  the  abscissa,  e,  is  the  strain  and 
the  ordinate,  a^W*,  Is  the  stress.  When  the  bar  is  free  of  body  forces, 
equilibrium  configurations  are  states  of  uniform,  l.e.  spatially  constant, 
stress  throughout  the  bar.  The  non-monotonlclty  of  the  stress-strain  re¬ 
lation  leads  to  an  Interval  (o  ,a_)  of  possible  values  of  equilibrium 

ct  p 

stress  for  which,  at  any  location  In  the  bar,  there  are  three  possible 
values  of  strain.  Note  that  In  the  Interval  (a  ,a.)  there  occurs  a 

Ot  p 

unique  stress  level,  called  the  Maxwell  stress  and  labeled  In  Fig.  1, 
with  the  geometric  property  that  the  two  closed  loops  in  the  stress-strain 
curve  determined  by  the  Maxwell  stress  have  equal  area.  In  using  a  tradi¬ 
tional  equilibrium  criterion  according  to  which  the  total  energy  In  the 
bar  Is  an  extremum,  we  distinguish  between  the  following  equilibrium  con¬ 
figurations.  An  equilibrium  configuration  Is  termed  stable,  neutrally 
stable,  or  unstable  with  respect  to  some  other  configuration  satisfying 
the  loading  conditions.  If  the  total  energy  of  the  equilibrium  configura¬ 
tion  Is  respectively  less  than,  equal  to,  or  greater  than  the  total  energy 
of  the  competing  configuration.  If  an  equilibrium  configuration  Is  stable 
with  respect  to  all  such  competing  configurations  we  say  It  is  absolutely 
stable;  If  It  Is  stable  or  neutrally  stable  with  respect  to  all  such 
competing  configurations,  we  say  it  Is  stable.  We  call  the  equilibrium 
configuration  metastable  If  It  is  not  stable,  but  If  all  competing 


configurations  with  respect  to  which  the  equilibrium  configuration  Is  unstabl 
have  strain  fields  that  are  not  uniformly  close  to  the  strain  field  of  the 
equilibrium  configuration.  Hence  a  metastable  equilibrium  configuration 
is  either  stable  or  neutrally  stable  with  respect  to  all  configurations 
whose  strain  fields  differ  only  Infinitesimally  from  the  strain  field  of 
the  equilibrium  configuration  everywhere  In  the  bar.  Finally  If  a  con¬ 
figuration  Is  neither  stable  nor  metastable  we  say  It  Is  unstable.  Con¬ 
sequently  the  type  of  stability  (stable,  metastable,  unstable)  depends  on 
the  size  of  the  class  of  configurations  with  which  the  equilibrium  con¬ 
figuration  successfully  competes  as  a  mlnlmlzer  of  the  total  energy. 
Stable  configurations  are  strong  relative  mlnlmlzers  of  the  energy  func¬ 
tional  In  that  competitors  need  only  satisfy  the  loading  condition,  so 
that  any  configuration  Is  an  eligible  competitor  provided  It  Is  In  the 
function  space  In  which  the  solution  Is  sought.  On  the  other  hand, 
metastable  configurations  are  merely  weak  relative  mlnlmlzers  of  the 
energy  functional,  since  the  number  of  competitors  is  reduced  by  requiring 
In  addition  that  the  strain  fields  In  the  competitors  be  uniformly  close 
to  the  strain  field  of  the  mlnlmlzer.  Let  <Y2<'f3  the  three 
Maxwell  strains  as  depicted  In  Fig.  1.  For  a  given  equilibrium  stress 
level  a.  It  is  shown  In  [1]  that  If  a<c^  (a>a^)  the  unique  config¬ 
uration  with  £<y.|  (e>Yj)  everywhere  in  the  bar  Is  absolutely  stable. 

If  a  lies  In  the  Interval  Caa»ag]  and  a<aY(a>ay)  there  is  one 
possible  strain  Ei  < Y-j ( E]  > y3)  and  another  e2  with  <*2ie2<Y3 
( y-j  sc?<6^).  It  Is  shown  that  a  configuration  of  the  bar  whose  strain 
field  Is  piecewise  uniform  and  takes  only  the  two  values  of  and 

Is  metastable.  If  a  lies  In  the  Interval  (a  ,oQ)  and  a<a  (a>a  ), 

a  p  Y  y 


-5- 


there  Is  a  third  possible  strain  e3  with  8-j  <  £3  <  otgj  any  configuration 
whose  strain  field  takes  the  value  e3  at  any  point  Is  unstable.  If 
0*0^,  there  are  no  absolutely  stable  configurations,  but  any  configura¬ 
tion  with  piecewise  uniform  strain  taking  only  the  values  and  y3  *s 
stable.  If  0*0^  and  c-yg  anywhere  In  the  bar,  the  configuration  Is 
unstable.  A  particle  P  of  the  bar  which  was  originally  at  x  will 
be  said  to  be  In  the. low-strain  (or  first)  phase  If  the  strain  e(x) 
corresponds  to  a  point  lying  below  the  first  ascending  branch  of  the 
stress-strain  curve  (Fig.  1),  so  that  Os£(x)sB^.  If  B^<e(x)<a2  . 

P  Is  said  to  be  In  the  unstable  phase,  while  P  Is  In  the  high-strain 
(or  second)  phase  If  e(x) 

Consider  a  bar  with  one  end  fixed  while  the  other  end  Is  loaded 
monotonl cal 1 y  and  quaslstatlcally.  The  resulting  equilibrium  states 
will  Involve  uniform  stress  at  each  Instant.  As  long  as  the  load  level 
Is  below  the  Maxwell  stress  one  would  expect  the  bar  to  equilibrate  in 
the  absolutely  stable  first  phase  state  with  a  uniform  strain  lying 
below  the  first  ascending  branch.  However,  once  the  load  exceeds  the 
Maxwell  stress  o^,  the  configuration  In  which  all  particles  are  In  the 
low-strain  phase  Is  only  metastable,  and  so  one  would  expect  to  eventually 
observe  strains  at  the  stabler  value  lying  below  the  second  ascending 
branch.  Transitional  equilibrium  configurations  would  then  be  character¬ 
ized  by  a  partitioning  of  the  bar  Into  a  number  of  co-existent  phases  of 
alternately  higher  and  lower  strains.  As  the  load  Is  Increased  still 
further,  the  regions  of  high  strain  would  grow  at  the  expense  of  the 
regions  of  low  strain,  leading  eventually  to  the  absolutely  stable  state 
In  which  the  bar  Is  everywhere  In  the  high- strain  phase.  Finally,  if 


one  then  conducts  an  analogous  program  of  unloading,  one  would  expect  to 
observe  portions  of  the  bar  In  the  low-strain  phase  at  some  load  below 
the  Maxwell  stress.  Hence  one  anticipates  a  hysteresis  loop.  The  time 
at  which  these  phase  transitions  occur  would  depend  on  the  actual 
physical  experiment  where  a  certain  unavoidable  noise  or  disturbance 
level  prevents  achieving  perfectly  smooth  loadings  and  unloadings.  If 
the  disturbance  level  is  high,  one  expects  that  the  bar  would  everywhere 
jump  to  the  stabler  phase  at_  the  Maxwell  stress.  However  if  it  is  low, 
one  expects  to  be  able  to  preserve  the  metastable  one-phase  states  almost 
until  that  phase  will  no  longer  sustain  the  given  equilibrium  stress. 
Consequently  the  size  of  a  measured  hysteresis  loop  would  be  an  indication 
of  the  disturbance  level  in  the  experiment. 

As  yet  we  have  not  excluded  the  possibility  that  different  locations 

in  the  bar  could  change  phase  at  different  equilibrium  stress  levels. 

For  a  homogeneous  bar  In  the  envisioned  loading  experiment,  one  would 

assume  that  all  transitions  from  low  to  high  strain  would  occur  at  the  same 

stress.  Hence,  In  equll  Ibrlum,  co-existence  phases  would  occur  only  at  some 

particular  value  of  stress  between  o  and  aQ.  The  choice  of  such  a  tran- 

o.  p  - 

sition  stress  under  loading  would  then  be  a  statement  about  the  disturbance 
level  of  the  system  and  amounts  to  choosing  a  particular  Inverse, 
e  *  (W*  )“”'  (a) ,  to  the  stress  response  o  =  W#(e).  Ericksen  [1]  showed  that 
when  this  transition  occurs  at  the  Maxwell  stress,  an  arbitrary  number  and 
placement  of  alternating  regions  of  hlgh-and  low-strain  phases  at  this 
stress  Is  consistent  with  either  a  specified  displacement  (hard)  end 
condition  or  a  specified  load  (soft)  end  condition.  However,  for  the 


hard  condition  ,  the  total  extent  of  each  phase  present  Is  uniquely 


determined  from  the  condition  which  specifies  the  total  length  of  the 
deformed  bar. 

An  analogous  situation  occurs  In  the  bending  of  an  Initially  straight, 
homogeneous,  Inextenslble  elastlr.a  with  a  non-convex  bending  energy 
density  [2],  In  order  to  facilitate  a  discussion  of  the  similarities  In 
the  bar  and  elastlca  problems,  we  shall  Introduce  nomenclature  from 
classical  Gibbsian  equilibrium  thermodynamics.  We  identify  a  conjugate 
pair  of  dependent  field  variables  for  each  system  under  consideration  In  such 
a  way  that  the  variable  which  appears  In  the  argument  of  the  energy 
density  will  be  called  the  extensive  variable,  while  the  first  derivative 
of  the  energy  density  with  respect  to  this  variable  will  be  called  the 
intensive  variable  [3],  For  the  bar  and  elastica  problems,  the  extensive 
variables  are  respectively  the  strain  and  the  curvature,  while  the  inten¬ 
sive  variables  are  respectively  the  stress  and  the  moment.  For  energy 
densities  of  the  type  studied  by  Erlcksen  [1],  Fig.  1  applies,  provided 
the  abscissa  Is  Identified  with  the  extensive  variable  and  the  ordinate 
Is  Identified  with  the  Intensive  variable.  Again  we  say  a  location  in 
the  body  Is  in  one  of  three  distinct  phases  depending  on  whether  the 
extensive  variable  Is  in  the  Interval  lying  below  the  first  ascending, 
descending,  or  second  ascending  branch  of  Fig.  1.  In  the  body,  the  ex¬ 
tensive  variable  may  be  discontinuous  across  certain  shock  surfaces . 

When  the  discontinuity  separates  values  of  the  extensive  variable  which 
lie  below  distinct  branches  of  the  curve  of  Fig.  1,  such  a  shock  surface 
Is  simultaneously  a  phase  boundary.  In  all  cases,  global  equilibrium 
requires  the  continuity  of  the  Intensive  variable. 


-8- 


Fosdlck  and  James  [2]  examine  two  problems  for  an  elastica  for  which 
Fig.  1  now  describes  the  moment-curvature  relation.  Equilibrated  bending 
moments  are  prescribed  in  one  problem,  while  the  other  is  one  of  pre¬ 
scribed  slope  difference  between  the  ends  of  the  elastica.  Each  problem 
is  one  of  "pure  bending"  in  that  equilibrium  configurations  are  those 
with  uniform  moment  throughout  the  elastica.  By  minimizing  the  total 
energy  subject  to  the  given  end  conditions,  we  distinguish  the  same  sta¬ 
bility  types  introduced  previously.^  Once  again  if  the  equilibrium  value 
of  the  intensive  variable  — in  this  case  the  moment— is  not  equal  to  o^, 
that  unique  equilibrium  configuration  in  which  the  extensive  variable - 
in  this  case  the  curvature -is  either  everywhere  less  than  y-j  or 
everywhere  greater  than  y3  is  absolutely  stable.  If  the  intensive 
variable  is  equal  to  a^,  there  are  no  absolutely  stable  configurations; 
but  equilibrium  configurations  in  which  the  extensive  variable  is  limited 
to  either  y-j  or  y3  are  stable.  This  case  includes  co-existent  phases 
in  which  the  deformed  elastica  consists  of  smoothly  connected  circular 
arcs  whose  curvatures  alternate  between  y-j  and  y3-  When  the  intensive 
variable  is  not  equal  to  c^,  such  coexistent  phases  between  the  two 
ascending  branches  of  Fig.  1  are  only  metastable.  In  every  case,  if  at  any 
point  in  the  body  the  value  of  the  extensive  variable  1  ies  below  the  descending 
branch,  the  system  is  unstable.  For  stable  equilibrium  configurations, 
co-existent  phases  occur  only  at  the  Maxwell  moment  a and  the  number 
of  alternating  regions  of  different  phase,  as  well  as  their  location,  is 
arbitrary.  Once  again,  for  the  hard  boundary  condition  -  prescribed 

Vosdick  and  James  in  [2]  use  a  different  terminology.  Our  stable  and 
metastable  equil ibrium  configuration  correspond  respectively  to  their 
Eulerian  and  weak  Eulerian  states. 
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termlnal  slope  difference  —  the  total  extent  of  the  elastica  in  each  phase 
Is  fixed.  No  such  restriction  applies  to  the  problem  with  the  soft  bound¬ 
ary  condition. 

For  the  problems  reviewed  above,  equilibrium  demands  that  the  in¬ 
tensive  variable  be  spatially  constant.  This  leads  to  highly  non-unique 
stable  equilibrium  configurations  when  that  constant  is  the  Maxwell 
value  o^.  In  more  complicated  problems,  the  intensive  variable  need 
not  be  independent  of  position.’  This  eliminates  much  of  the  arbitrariness 
In  the  stable  equilibrium  configurations  with  co-existent  phases.  For 
example,  bars  where  each  particle  has  the  same  type  of  stress  response  as 
in  Fig.  1.,  but  with  material  and  geometrical  inhomogeneities,  admit  the 
possibility  that  the  transition  stress  may  vary  throughout  the  bar.  In 
addition  if  body  forces  are  present,  equilibrium  configurations  need  not 
be  configurations  of  spatially  constant  stress.  Consequently  the  stress 
in  the  bar  will  coincide  with  the  transition  stress  only  at  certain  loca¬ 
tions.  As  shown  by  James  [4],  this  permits  certain  conclusions  to  be 
drawn  about  the  number  and  location  of  equilibrium  phase  boundaries. 

Another  example  is  provided  by  an  elastica  with  one  end  fixed  and 
subject  to  an  axial  compressive  force  at  the  other  end.  Here  buckling 
becomes  a  further  potential  source  of  non-uniqueness.  For  nonlinear 
moment-curvature  relations,  the  number  of  available  stable  buckled 
configurations  will  in  general  increase  at  certain  values  of  the  applied 
load.  For  each  such  buckled  state,  the  associated  finite  deformation 
provides  a  variable  moment-arm  for  the  applied  load,  so  that  equilibrium 
configurations  will  not  be  those  with  a  uniform  moment  field.  James  [5] 
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ha  s  studied  buckled  configurations  when  the  moment -curvature  relation 
Is  as  In  Fig.  1.  He  has  shown  that  when  the  applied  load  Is  sufficiently 
large,  candidates  for  stable  "first  mode"  buckled  configurations  Involve 
a  region  of  high-curvature  phase  at  the  fixed  end  of  the  elastlca,  while 
the  remaining  portion  Is  In  the  low-curvature  phase.  The  location  of 
the  phase  boundary  Is  precisely  determined,  being  further  from  the  fixed 
end  for  larger  applied  loads. 

A  three-dimensional  example  Is  a  problem  of  the  twisting  of  an  In¬ 
compressible  homogeneous  elastic  tube  with  a non-convex  energy  density,  as 
studied  by Abeyaratne  [6].  By  seeking  radially  symmetric  solutions, 
Abeyaratne  reduces  the  problem  to  a  one-dimensional  one.  Here  the  In¬ 
tensive  variable  Is  the  shear  stress,  which  Is  not  In  general  constant  for 
equilibrium  configurations.  It  is  found  that  minimizers  of  the  total  stored 
energy1 correspond  to  a  unique  configuration  for  all  applied  twisting 
angles.  Some  of  these  solutions  are  smooth,  while  others  Include  a 
single  circular  phase  boundary. 

Knowles  and  Sternberg  [7]  have  examined  the  relationship  between 
the  loss  of  elllptlclty  of  the  governing  equations  and  the  emergence  of 
equilibrium  shocks  In  the  two-dimensional  problem  of  plane  finite  elasto- 
statlcs.  It  Is  shown  that  a  necessary  condition  for  the  existence  of  a 
straight  shock  separating  two  homogeneous  plane  deformations  Is  that  the 
equilibrium  equations  must  suffer  a  loss  of  strong  elliptlclty  at  some 

1 

Abeyaratne  restricts  competitors  to  be  symmetric  and  also  places  an  L2 
norm  restriction  on  the  closeness  of  the  shear  between  minimlzer  and  com¬ 
petitors.  This  restriction  is  weak  enough  not  to  unduly  limit  the  class 
of  competitors  and,  as  in  [2,4,5]  where  there  are  no  norm  restrictions, 
the  jump  in  the  Intensive  variable  must  occur  at  the  Maxwell  value. 
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homogeneous  deformation  which  Is  a  convex  linear  combination  of  the  two 
plane  homogeneous  deformations  separated  by  the  shock.  This  fits  In 
nicely  with  the  one-dimensional  problems  reviewed  above,  where  It  was 
found  that  In  stable  and  metastable  equilibrium  configurations  the  values 
of  the  extensive  variable  on  opposite  sides  of  a  phase  boundary  are  sepa¬ 
rated  by  an  unstable  branch  of  the  constitutive  relation. 

The  convexity  of  the  relevant  energy  density  plays  a  central 
role  with  respect  to  the  question  of  the  existence  of  equilibrium 
solutions  which  minimize  the  energy  functional.  For  the  elastic  bar. 

Ball  [8]  demonstrates  that  If  all  specified  displacement  problems  are 
to  have  continuously  differentiable  minlmlzers  of  the  energy  functional, 
then  the  corresponding  strain  energy  density  must  necessarily  be  convex. 
When  It  Is  not,  certain  problems  of  specified  displacement  will  have 
minlmlzers  only  If  they  admit  discontinuous  strains.  Conversely,  con¬ 
vexity  of  the  strain  energy  Is  also  a  sufficient  condition  for  the 
existence  of  minlmlzers. 

Gurtln  and  Temam  [9]  have  considered  the  minimization  problem  for 
finite  elastostatic  anti-plane  shear.  Here  the  extensive  variable— the 
shearing  stress -and  the  Intensive  variable -the  shearing  strain -are 
related  as  in  Fig.  1.  They  consider  an  auxiliary  minimization  problem 
governed  by  a  strain  energy  density  which  is  the  lower  convex  envelope 
of  the  original  strain  energy  density.  Then  the  "convexified  problem"  Is 
governed  by  a  relation  between  shearing  stress  and  shearing  strain  In 
which  the  Maxwell  line  replaces  that  portion  of  the  original  curve  which 
lies  above  the  metastable  and  unstable  Intervals  of  shearing  strain. 
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Yl<e<Y3‘  ^n^e  the  original  problem,  the  convexlfled  problem  Is 
guaranteed  to  have  a  mlnlmlzer.  The  minimum  value  of  the  convexlfled 
energy  functional  Is  the  Inflmum  of  the  energy  functional  of  the  original 
problem  evaluated  at  all  configurations  which  satisfy  the  boundary  con¬ 
ditions.  Conversely  If  mlnlmlzers  of  the  original  problem  exist,  they 
will  be  mlnlmlzers  of  the  convexlfled  problem. 

For  three-dimensional  finite  el astostatlcs,  Abeyaratne  [10]  has  shown  that 
the  equilibrium  shock  conditions  are  given  naturally  by  the  Welerstrass- 
Erdmann  corner  conditions  when  piecewise  smooth  functions  are  admitted 
Into  contention  as  possible  mlnlmlzers  of  the  energy  functional.  The 
traction  continuity  condition  arises  from  considering  variations  of  the 
displacement  components.  If  one  then  considers  variations  of  the  Indepen¬ 
dent  variables,  the  spatial  coordinates,  another  vectorial  shock  condition 
Is  obtained.  Continuity  of  the  tangential  displacement  derivatives  across 
the  shock  guarantees  that  all  but  one  component  of  this  shock  condition 
Is  satisfied  automatically.  One  Is  then  left  with  a  supplemental  shock 
condition  requiring  an  energy-like  quantity  to  be  continuous  across  the 
discontinuity  surface.  For  the  equilibrium  bar  problem  of  Erlcksen  [1], 
this  condition  Is  satisfied  for  a  material  with  shocks  If  and  only  If 
the  constant  equilibrium  stress  In  the  bar  l£  the  Maxwell  stress. 

Knowles  and  Sternberg  [11,  12]  and  Abeyaratne  [13]  have  examined 
certain  anti-plane  shear  crack  problems  with  a  non-convex  energy  density 
and  have  found  equilibrium  solutions  with  shocks  Issuing  from  the  crack 
tips  [11,12]  or  from  points  on  the  crack  surface  [13]  and  terminating 
In  the  Interior  of  the  body.  These  solutions  do  not  satisfy  the  sup¬ 
plemental  shock  condition.  If,  however,  one  considers  a  one-parameter 
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family  of  equilibrium  states  corresponding  to  a  quasistatic  loading,  it 
can  be  shown  that  across  the  shocks  appearing  in  the  corresponding 
quasistatic  solutions  of  [11,12,13]  the  jump  in  the  supplemental  quantity 
is  of  a  particular  sign  (with  respect  to  the  direction  with  which  the 
quasistatic  shock  is  moving).  Thus  although  the  supplemental  shock 
condition  Is  not  satisfied,  a  supplemental  shock  inequality  is  maintained. 
Moreover  this  Inequality  assures  the  dissipative  character  of  the  shock  in 
the  sense  that  It  precludes  quasistatic  shocks  which  allow  the  body  to 
store  elastic  energy  faster  than  the  rate  at  which  work  is  being  done. 

This  disslpatlvity  inequality  was  originally  proposed  as  an  admissibility 
criterion  for  elastic  shocks  by  Knowles  and  Sternberg  [7]  and  Knowles  [14]. 
It  is  analogous  to  the  entropy  Inequality  for  gas-dynamical  shocks. 

For  the  quasi  statically  loaded  bar,  this  disslpatlvity  inequality 
allows  the  emergence  of  shocks  only  If  the  stress  equals  or  exceeds  the 
Maxwell  stress.  Similarly  Posdick  and  James  [2]  find  that  the  metastable 
equilibrium  configurations  for  the  pure  bending  problem  of  the  elastica 
satisfy  an  analogue  of  the  dissipation  inequality  across  the  curvature 
discontinuities.  It  thus  appears  that  metastable  quasistatic  proces¬ 
ses  are  associated  with  dissipative  and  hysteretlc  phenomena,  while 
absolutely  stable  quasistatic  processes  are  dissipation-free. 

Of  course,  the  quasistatic  processes  envisioned  above  are  not  in 
reality  dynamic  processes  at  all,  but  rather  one-parameter  families  of 
equilibrium  states.  Although  it  Is  logical  to  identify  the  history 
parameter  with  time,  such  an  approach  neglects  the  effect  of  inertia  and 
the  consequent  wave  propagation  phenomena  associated  with  dynamical 
systems.  Studying  the  emergence  and  evolution  of  different  phases  in  a 


fully  dynamic  theory  would  seem  to  be  prudent  both  for  experimental 
reasons  and  to  better  understand  the  approach  to  an  equilibrium  consisting 
of  co-existent  phases.  The  latter  would  presumably  help  to  resolve 
some  of  the  questions  unanswered  by  the  static  theory.  In  particular 
the  ultimate  location  of  equilibrium  phase  boundaries. 

The  following  brief  discussion  of  the  dynamic  behavior  of  systems 
with  non-convex  equilibrium  energy  densities  will  be  limited  to  one¬ 
dimensional  problems  for  bars.  Even  then  there  Is  a  wide  variety  of 
dynamic  problems  which  may  be  related  to  any  one  particular  quasi  static 
problem,  since  non-equilibrium  processes  can  be  governed  by  an  infinity 
of  constitutive  assumptions,  all  of  which  reduce  in  an  equilibrium  setting 
to  the  same  purely  elastic  constitutive  relation. 

Determining  the  ultimate  location  of  phase  boundaries  makes  visco¬ 
elastic  theories  particularly  appropriate,  since  the  Inherent  dissipation 
of  a  viscoelastic  dynamical  process  can  lead  to  an  asymptotic  static  state 
for  large  time.  Dafermos  [15]  has  studied  an  initial-boundary  value 
problem  for  a  special  class  of  viscoelastic  bars  for  which  the  constitu¬ 
tive  law  in  the  static  case  Is  consistent  with  a  non-monotonic  equil ibri urn 
stress-strain  relation.  He  considers  an  initially  deformed  viscoelastic 
bar  which  Is  released  with  some  initial  velocity  and  is  subsequently  free 
from  body  forces  and  end  loads.  It  Is  shown  that  if  the  viscosity  is 
bounded  above  zero  the  bar  will  asymptotlcaT'y  approach  some  stress-free 
equilibrium  state.  Especially  significant  In  light  of  our  review  of  the 
static  theory  is  that  the  asymptotic  deformation  gradient  will  in  general 
be  discontinuous,  possibly  even  unbounded.  Moreover,  It  Is  not  completely 
apparent  which  of  the  various  possible  equilibrium  configurations  is 


ultimately  approached. 

The  motion  of  a  viscoelastic  bar  Is  typically  governed  by  a  partial 
differential  equation  of  order  three  or  more.  In  which  case  the  well- 
developed  mathematical  techniques  for  second-order  partial  differential 
equations  are  unavailable.  However  the  dynamics  of  a  purely  elastic  bar 
are  governed  by  a  second-order  quaslllnear  partial  differential  equation 
of  mixed  type,  being  hyperbolic  when  the  strain  lies  below  an  ascending 
branch  of  the  stress-strain  curve,  and  elliptic  when  the  strain  lies 
below  a  descending  branch.  As  naturally  occurring  Initial  and  boundary 
conditions  fordynamlcal  problems  lead  to  well-posed  problems  for  hyper¬ 
bolic  partial  differential  equations  but  111  posed  ones  for  elliptic 
equations,  the  dynamic  theory  Itself  provides  Impetus  to  seek  solutions 
of  the  dynamic  problem  which  have  strains  confined  to  the  statically 
stable  and  metastable  phases.  With  such  a  restriction  the  governing 
equation  Is  equivalent  to  a  hyperbolic  system  of  two  first-order  ordinary 
differential  equations  and  so  the  theory  of  Rlemann  Invariants  and 
characteristic  curves  In  the  x-t  plane  is  at  one's  disposal.  Space- 
time  curves  across  which  the  strain  jumps  from  one  branch  to  the  other 
are  the  propagating  phase  boundaries. 

To  appreciate  the  specific  dynamical  phenomena  associated  with 
non-monotonlc  stress-strain  relations  It  is  helpful  to  first  consider  the 
simpler  theory  for  the  monotonic  case.  A  good  account  of  this  theory  may  be 
found  In  Courant  and  Friedrichs  [16].  It  Is  well  known  [17]  that  solu¬ 
tions  for  hyperbolic  systems  with  arbitrarily  smooth  Initial  data  are 
not  guaranteed  to  be  globally  smooth  for  all  time;  the  elastic  bar  Is 
no  exception.  The  type  of  singularities  encountered  In  this  setting 


are  associated  with  the  Intersection  of  characteristic  curves  at  some 
finite  distance  from  the  curve  of  Initial  data.  At  such  points  the 
second  derivatives  of  the  displacement  become  unbounded.  The  remedy  Is 
to  seek  weak  solutions  to  the  governing  equations  by  Introducing  shock 
curves  across  which  the  first  derivatives  of  the  displacement  are  dis¬ 
continuous,  but  their  jumps  satisfy  two  shock  conditions  which  are 
necessary  for  displacement  continuity  and  momentum  balance.  As  there 
may  be  many  such  weak  solutions,  an  "entropy  condition",  due  to  Lax  [18], 
may  be  Introduced  to  select  solutions  of  physical  significance.  This 
has  the  effect  of  confining  the  shock  speed  to  the  Interval  between  the 
different  acoustic  speeds  on  each  side  of  the  shock. 

It  must  be  emphasized  that  these  types  of  shocks  are  to  be  distin¬ 
guished  from  phase  boundaries,  since  the  latter  are  not  even  admitted 
Into  a  bar  theory  with  a  monotone  stress-strain  relation.  For  the  non¬ 
monotonic  theory  both  types  of  shocks  become  possible.  By  a  conventional 
shock  we  shall  mean  a  shock  which  arises  from  the  Intersection  of  character 
Istlc  curves  associated  with  strains  lying  below  a  single  branch  of  the 
stress-strain  relation.  Shock  curves  which  separate  distinct  phases  are 
referred  to  as  phase  boundaries.  The  term  "shock"  without  an  obvious 
referent  will  be  used  In  discussing  features  common  to  both  types.  For 
example, as  weak  solutions  require  the  same  relations  to  hold  across  all 
curves  with  discontinuous  strain,  the  same  "shock"  conditions  apply  in 
both  cases. 

The  central  role  of  monotonicity  and  convexity  conditions  in  general 
hyperbolic  systems  of  equations  has  been  explored  by  Lax  [18],  By 
applying  the  methods  used  In  [18]  to  the  equations  for  an  elastic  bar 
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with  a  stress-strain  relation  as  In  Fig.  1,  James  [19]  has  demonstrated 
the  local  existence  of  classically  smooth  fields  about  an  assigned 
moving  phase  boundary.  James  also  treats  two  specific  Initial  value 
problems.  In  the  first, a  value  of  strain  Is  chosen  from  each  of  the  two 
Intervals  lying  below  ascending  branches  of  Fig.  1.  One  value  Is  specified 
along  the  positive  x  axis  and  the  other  along  the  negative  x  axis.  Hence 
the  bar  Initially  contains  two  phases.  Two  appropriately  selected  ve¬ 
locities  complete  the  specified  initial  data.  For  this  problem,  two 
one-paramteter  families  of  solutions  are  found,  each  of  which  contains 
a  single  constant-velocity  phase  boundary  emerging  from  the  point  of 
initial .  strain  discontinuity.  The  second  Initial  value  problem  treated 
In  [19]  Involves  a  constant  value  of  strain  associated  with  one  of  the 
ascending  branches,  and  zero  velocity  specified  Initially  all  along  the 
x  axis.  In  addition  to  the  obvious  static  solution  In  which  the  Initial 
stable  or  metastable  state  persists  for  all  time,  James  establishes  a 
two-parameter  family  of  solutions  which  contain  both  a  constant  velocity 
phase  boundary  moving  to  the  right  and  to  the  left.  The  region  between 
the  phase  boundaries  experiences  a  constant  strain  lying  below  the  other 
ascending  branch  of  the  stress-strain  curve.  James  points  out  that  this 
second  problem  Is  naturally  associated  with  necking  phenomena, in  that 
phase  boundaries  may  spontaneously  arise  In  the  Interior  of  the  bar.  It 
would  however  seem  that  achieving  the  initial  condition  of  spatially 
constant  non-zero  strain  would  Itself  require  a  quasistatic  process.  Also 
the  constant  speed  of  the  two  phase  boundaries  depends  on  the  doubly 
Infinite  spatial  domain,  so  that  subsequent  boundary  effects -such  as 
those  due  to  a  loading  device  or  due  to  waves  reflected  from  the  ends  of 
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the  bar— are  not  taken  Into  account. 

James  also  addressed  the  question  of  admissibility  of  solutions  to 
the  dynamic  problem  by  alternately  considering  static,  viscous  and 
entropy  rate  criteria.  Particularly  Interesting  Is  the  "viscosity"  cri¬ 
terion,  according  to  which  an  Inequality  to  be  satisfied  across  shocks  Is 
deduced  by  considering  solutions  which  are  limits  of  viscoelastic  solutions 
as  the  viscocity  vanishes.  Fora  phase  boundary  In  which  the  traction  discon¬ 
tinuity  Is  small— so  that  the  phase  boundary  Is  moving  slowly  and  sepa¬ 
rates  a  region  of  statically  stable  phase  from' a  region  of  statically 
metastable  phase— the  viscoelastic  criterion  requires  that  the  phase 
boundary  shall  move  In  the  direction  that  converts  statically  metastable 
phase  to  statically  stable  phase.  Hence  when  the  traction  discontinuity 
Is  small,  both  the  viscoelastic  criterion  of  James  and  the  dissipation 
criterion  of  Knowles  allow  and  exclude  the  same  phenomena.  Indeed  as 
the  traction  discontinuity  vanishes,  the  phase  boundary  ceases  Its  motion 
and  the  distinction  between  the  two  criteria  vanishes.^  For  the  first 
Initial  value  problem  considered  by  James,  the  viscoelastic  criterion 
allows  only  one  of  the  two  families  of  solutions  found.  For  the  second 
problem,  the  occurrence  of  the  double  phase  boundaries  Is  possible  only 
If  the  constant  Initial  strain  Is  a  statically  metastable  solution. 


^Compare  (5.39)  of  [19]  and  (4.19)  of  [14]. 
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2.  Governing  Equations  for  the  Dynamically  Loaded  Elastic  Bar. 
Characteristic  Theory  when  a  Single  Phase  is  PresenTT 


A.  Formulation  of  the  Problem 

We  consider  a  homogeneous,  semi-infinite  elastic  bar  which  In  the 
reference  configuration  occupies  the  Interval  -  m<xsO.  Motions  are 
described  by  a  mapping 

y(x,t)  -x  +  u(x,t)  ,  (2.1) 


where  y(x,t)  Is  the  coordinate  at  time  t  of  a  particle  which  Is  at  x 
In  the  reference  configuration;  u(x,t)  Is  the  displacement.  We  assume 
the  reference  configuration  to  be  undeformed. 

In  the  fully  three  dimensional  theory  of  elasticity  the  strain  energy 
density  W  is  a  function  of  the  deformation  gradient  tensor  £  .  When 
the  above  deformation  Is  expressed  In  a  Cartesian  reference  frame  J  with 
the  unit  vector  ^  aligned  with  the  bar,  £  Is  described  by  means  of  its 
matrl x  ; 


1  +  u'  0  0 

0  1  0 

0  0  1 


(2.2) 


A  prime  will  In  general  denote  differentiation  with  respect  to  the  argu 
ment  and  specifically  with  respect  to  x  If  there  are  multiple  Indepen 
dent  variables.  The  Plola  stress  tensor  ^  is  given  by 

Z  *  ^  * 


(2.3) 
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For  the  longitudinal  deformation  (2.1)  of  an  elastic  bar,  the  strain 
energy  density  becomes  a  function  of  the  displacement  gradient -or 
strain— which  we  denote  by 

e(x,t)  au'(x,t),  (2.4) 

so  that 

W(JF)  -  W(e) .  (2.5) 

It  Is  assumed  that  for  this  type  of  deformation  W(e)  Is  such  that  the 
only  nonzero  stress  Is  the  normal  component  along  the  axis  of  the  bar, 
which  shall  be  denoted  by  a  and  according  to  (2.3)  Is  given  by 

a(e)  *  Vf' (e)  .  (2.6) 

The  longitudinal  velocity  In  the  bar  Is 

v(x,t)  -  y(x,t)  »  u(x,t)  ,  (2.7) 

where  the  superposed  dot  denotes  time  dlfferentatlon. 

With  no  body  forces  the  only  nontrivial  equation  of  motion  Is  that 
associated  with  the  balance  of  momentum  In  the  x  direction: 

pU  ■  W* (u' )  »  W*  (u' )u*  ,  (2.8) 

where  p,  a  positive  constant.  Is  the  density  of  the  undeformed  bar. 

For  a  bar  that  Is  Initially  undisturbed  and  at  rest. 


u(x,0)  =  0,  v(x,0)  =  0, 


-  0D  <XiO  . 


(2.9) 
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We  shall  consider  two  different  types  of  loading  device.  In  the 
hard  device,  the  bar  is  subject  to  a  controlled  end  displacement  uQ(t), 
so  that 


u(Q,t)  *  uQ(t)  ,  tiO  .  (2.10) 

The  other,  a  soft  device,  subjects  the  bar  to  an  end  load  aQ ( t ) .  Here 
the  appropriate  boundary  condition  is 


W'(u'(0,t))  »  oQ(t)  ,  t*0.,  (2.11) 

We  shall  seek  solutions  u(x,t)  which  are  twice  continuously  dif¬ 
ferentiable  for  x^O,  tiO  except  on  at  most  a  finite  number  of  smooth 
shock  curves,  x*s^(t).  Where  u  is  twice  continuously  differentiable, 
we  require  that  (2.8)  be  satisfied.  Across  a  shock  curve  —  x  =  s(t) ,  say- 
we  shall  require  that  the  displacement  shall  remain  continuous  and  that 
a  global  balance  of  momentum  be  maintained.  This  will  hold  provided  the 
following  shock  conditions -sometimes  called  Rankine-Hugoniot  conditions  — 
are  satisfied  [20]: 


s[u']  +  [u]  =  0  , 
ps[u  ]  +  [W*  (u' )]  =  0  , 


(2.12) 


where  [  ]  denotes  the  jump  across  the  shock  (e  g. 

[u']  =  u'  (s(t)t  t)  -  (/  (s(t)~,  t)). 


B.  Comments  on  the  Formulation 

Many  of  the  equilibrium  studies  cited  in  the  previous  section 
[1,2, A, 5, 10]  were  formulated  in  a  variational  setting.  It  should  be 


mentioned  that  this  problem  too  has  a  variational  formulation;  one 
seeks  extremals  to  the  appropriate  action  functional,  namely  the  time 
Integral  of  the  Lagrangian  of  the  system.  This  yields  (2.8)  as  the 
Euler-Lagrange  equation.  Moreover,  the  Weierstrass-Erdmann  corner  condi¬ 
tion  which  arises  from  varying  the  strain  and  velocity  lead  to  the  shock 
condition  which  expresses  the  balance  of  momentum  (2.12)2*  Variations 
with  respect  to  x  and  t  yield  two  more  corner  conditions  which  can 
be  reduced  to  (2.12)^,  and  an  extra  condition  which  requires  that  an 
energy-like  quantity  be  conserved  across  a  discontinuity  curve.  For 
static  solutions  this  extra  condition  reduces  to  the  no-dissipation  con¬ 
dition  derived  in  [10]. 

It  can  also  be  shown  that  (2.12)  are  the  discontinuity  conditions 
associated  with  the  canonical  weak  formulation  of  this  problem  [19]. 
Whitham  [21]  points  out  that  different  weak  formulations  may  yield  the 
same  partial  differential  equations  in  regions  where  the  fields  are  suf¬ 
ficiently  smooth,  but  may  lead  to  different  shock  conditions.  Realizing 
this,  it  is  not  difficult  to  discover  that  smooth  solutions  to 

(J.pG2  +  W(u<))  =  Uw'(u'))  (2-13 

are  solutions  to  (2.8)  where  u^O.  If  u  is  not  differentiable,  the 
canonical  weak  formulation  of  (2.13)  yields  (2.12)^  and  the  extra  energy¬ 
like  shock  condition,  but  does  not  yield  (2.12)2.  Hence  we  defer  to  the 
physics  of  the  problem -in  this  case  the  balance  of  momentum-to  require 
(2.12)  as  the  only  shock  conditions  which  solutions  must  satisfy. 


C.  Smooth  Solutions  with  One  Phase  Present 


We  shall  consider  materials  with  a  non-convex  strain  energy  density 
W  and.  In  particular,  we  assume  that: 


(1) 

(11) 

(111) 


and 


W  Is  twice  continuously  differentiable  on  [0,co  ), 
W'(O)  *  0,  W'(e)  >0  for  e  >0  , 


there  exists  8-j ,  a2  with  B-j  <a2  and 
W"  (8^  »  W"  (otg)  *  0  , 

W*  (e)>0  for  0ie<6-j  ,  e>a2  • 
W*  (e)  <  0  for  8-j  <  £  <  a2  ; 


( 1 v )  11m  W* (e)  «  ©  , 
e-*-ao 


(2.1< 


(2.1! 


(2.11 


whence  there  exists  82>a2  with 

W*  (82)  ■  W't^)  ■  Og  .  (2.1  j 

This  leads  to  a  stress-strain  relation  of  the  type  discussed  In  the  pre¬ 
ceding  section  and  depicted  In  Fig.  1. 

When  e  lies  In  the  Interval  0se<6.|  or  e>a2  ,  the  governing 
equation  (2.8)  Is  hyperbolic;  It  Is  elliptic  when  8-|  <e<a2.  Conse¬ 
quently  If  we  assume  that  the  loading  conditions  (2.10)  or  (2.11).  lead 
to  a  first  phase  strain  field  with  0se<6^  for  all  xsO  and  tuO, 
then  any  smooth  solution  to  (2.8),  (2.9)  and  (2.10)  or  (2.11)  can  be 
found  by  the  theory  of  Rlemann  Invariants  and  characteristic  curves.  A 
full  account  of  this  theory  may  be  found  In  [16].  Define  the  local 


'  o  •: 


acoustic  speed  by 


,  0ses81  , 


(2.18) 


and  let 


fc, 

<j»(e)  -  |  c(s)d 8  ,  0ses61 


(2.19) 


From  (2.18)  It  follows  that  $(e)  is  str+ctiy  increasing  on  0se<8^  . 
Let  4>(z)  be  the  Inverse  of  $  defined  on  0sz<<|)(81).  Then  by  (2.18) 
and  (2.19),  $  Is  differentiable  with 


✓  (z)  -  (lw*  (*(z)))-*  . 


(2.20) 


We  note  that  if  W  €CN([0,81 )) ,  then  c  €  CN"2([0,81 )),  <p  €  CN"1([0,81 )), 
and  $6^  1  ( [0,<|>( ) ) ) .  One  finds  that  (2.8)  Is  equivalent  to 


0  on 


c(e)  , 


(2.21) 


-  0 


on  3?  ’  _c(e) 


These  equations  in  turn  imply  the  existence  of  two  families  of  characteristic 
curves.  In  the  first  family  (2.21)^,  each  member-curve  has  positive  slope 
in  the  x-t  plane.  We  parametrize  this  family  by  the  variable  C+  . 

The  second  family  (2.21)2*  whose  member-curves  have  negative  slope,  we 
parameterize  by  C_  .  It  now  follows  from  integrating  (2.21)  that 
there  exist  functions  K+  and  K"  such  that 


*  »  *  i  «•  •' ■  •  *  •  •  •  %  %  *, 

J  V  V.V  --  V  mV  -*  /  •  _  /.  nm_  ,  '  , 


•  V*V*  .*•  /• 


(2.22) 


v-$(e)  a  K  (C+)  on  curves  of  constant  C+  given  by  ^  *  c(e), 

► 

v  +  $(e)a  K"(C_)  on  curves  of  constant  C_  given  by  &  «  -c(e). 

J 

Since  ^|->0  on  each  curve  C+,  all  curves  C+  originate  on  t«0, 
xsO  and  terminate  on  x«0,  taO.  On  the  other  hand  curves  C_  origi¬ 
nate  on  either  t*0,  xsO  or  x*0,  tiO  and  propagate  into  x<0, 
t  >  0.  Evaluating  K+  on  t«0,  xsO  one  infers  from  (2.4),  (2.9)  s 
(2.19)  and  (2.22)^  that  K+  =  0,  whence 

v  -  <fr(e)  .  (2.23) 

Then  on  each  curve  C_,  (2.22)  and  (2.23)  yield 

va^K'(Cj  ,  e  -  *(^K'(CJ)  .  (2.24) 

Moreover  (2.22)2  an<*  ^.24)g  show  that  each  curve  C_  is  a  straight  ray. 

If  we  further  assume  that  a  unique  ray  C_  passes  through  each  point 
(x,t),  xiO,  tsO,  then  there  is  a  simple  geometrical  construction  of 
the  solution.  Consider  first  the  soft  device  governed  by  (2.11).  The 
first  assumption  0<e<6j,  together  with  (2.15),  implies  that  the  end¬ 
load  induces  a  unique  end-strain  e(0,t)  3 eQ(t)  such  that 

W#(eQ(t))  *  aQ(t)  ,  tiO  ,  (2.25) 

and 

0seo(t)<61  .  (2.26) 

By  (2.15)  and  (2.26)  it  is  possible  to  solve  (2.25)  for  e  (t)  if  and 


0*oQ(t)  <*(6,)  «crB 


t  *0 . 


(2.27) 


9 


Consider  the  cay  C_  originating  at  x»0  and 
which  we  denote  by  C_(x),  has  slope  ^  •  -c(eq(t) ) . 
a  particular  (x,t)  lies  on  C_(t)  If  and  only  If 


t  ■  x .  This  ray. 
Consequently 


x  *  -c(£Q(x))(t  -  t).  (2.28) 

Since  e  and  v  are  constant  on  each  ray  C_,  It  follows  that  when  (2.28) 
holds,  e*e0(x).  Also  by  (2.23)  It  follows  that  v  *  $(e  (x)). 

Now  consider  a  particular  (x,t)  which  lies  on  a  ray  C_  originat¬ 
ing  on  t*0,  xsO,  so  that 

x  < -c(eo(0))t.  (2.29) 

The  constancy  of  e  and  v  on  C_  rays  together  with  the  Initial  condi¬ 
tion  (2.9)  then  Indicates  that  e*0  and  v»0.  Physically,  the  rays  C_ 
propagate  Information  about  the  end  loading  with  speed  c(e).  Consequently 
(2.29)  characterizes  the  reglonof  the  bar  In  front  of  the  leading  disturbance. 
Summarizing  these  results,  we  have: 


(1)  e  *  0 ,  v»0  for  x<-c(EQ(0))t 

(11)  £»E0(x(x,t)),  V  *  4>(  Eg(x(  X  »t )  )  for  Xi-c(Efl(0))t 


(2.30) 


where  x(x,t)  Is  given  Implicitly  by  (2.28).  We  shall  refer  to  (2.30) 
as  the  formal  solution  for  a  material  in  the  first  phase  loaded  by  a  soft 
device.  The  geometry  of  the  rays  C  is  depicted  In  Fig.  2. 
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Recall  that  In  arriving  at  (2.30)  we  made  two  assumptions,  the 
first  being  that  0se<8^  throughout  the  bar  for  all  time,  and  the 
second  that  a  unique  ray  passes  through  each  (x,t)  In  x<0,  t>0. 

We  now  examine  what  these  assumptions  entail.  According  to  (2.30),  the 
bar  will  be  everywhere  In  the  first  phase  with  e  <  8-j  If  the  same  Is 
true  of  the  Induced  end-strain  eQ(t).  Hence  (2.27)  is  both 
a  necessary  and  sufficient  condition  for  our  first  assumption.  We  note 
that  If  (2.27)  Is  not  satisfied,  then  by  (2.11)  and  (2.15)  the  end  load 
oQ(t)  cannot  be  sustained  by  an  end-strain  less  than  the  value  B-| . 
Consequently, the  formal  representation  (2.30)  would  then  Itself  be  mean¬ 
ingless  since  c(e)  and  <p(e)  are  only  defined  for  e^8-|  -  Indeed  an 
attempt  to  extend  the  definitions  (2.18)  and  (2.19)  Into  the  statically 
unstable  strain  Interval  8-j<e<<X2,  results  In  c(e)  and  <f>(e)  becoming 
Imaginary.  Most  of  this  paper  will  be  concerned  with  the  consequences  of 
smoothly  loading  the  bar  In  such  a  way  that  eventually  (2.27)  no  longer 
holds.  In  this  endeavor  It  will  be  convenient  to  have  (2.30)  available 
to  us  when  €<8^  consequently  It  Is  essential  that  we  also  examine  our 
second  assumption. 

Denote  the  disturbed  and  tranquil  regions  of  the  bar  by  Rq  and 
fty  where 


Rq  3  ((x,t)|t^0,  -c(eQ(0))t  *x  *0  j  , 

Rt  =  ((x,t)  |t  *0,  x  < -c(eo(0))t]  . 


(2.31) 


In  Ry,  all  the  characteristic  rays  C_  originate  on  x<0,  t=0,  so  that  by 
(2.9),  (2.4)  and  (2.21 )2  they  all  have  slope  ^  *  -c(0).  This  guarantees 
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that  the  rays  C_  form  an  unambiguous  cover  for  R^;  fn  other  words 
they  span  Rj-  and  do  not  intersect  each  other. 

In  Rq,  the  rays  C_  will  also  form  an  unambiguous  cover,  provided 
(2.28)  is  uniquely  Invertible  for  r(x,t)  whenever  (x,t)  €Rjj.  From  the 
construction  of  the  rays  C_  as  depicted  in  Fig.  2.,  this  will  occur 
provided  c ( ( t ) )  Is  a  continuous  and  monotonlcally  decreasing  function 
of  t  for  TiO.  If  oQ(t)  Is  continuous  or  differentiable,  It  follows 
from  (2.25),  (2.4)  and  (2.14)  that  eQ(t)  has  the  same  degree  of 
smoothness.  When  aQ(t)  is  differentiable 


cQ(t)  ■  o0(t)/lT  (eQ(t)) . 

Suppose  then  that 

oQ(t)  €  C^([0,®  ) )  , 
W(e)  €C3  ([0,6-j))  . 


(2.32) 


■\ 


(2.33) 


It  then  follows  that  c(£q(t))  Is  monotonlcally  decreasing  provided 


d 

TTF 


c(e _(t)) 


By  (2.15)  this  is  equivalent  to 


(2.34) 


W^(eo(t))oo(t)s0  ,  tiO  .  (2.35) 

Since  the  bar  Is  loaded  from  an  undeformed  rest  configuration,  it  is 
natural  to  consider 

aQ(t)  iO  ,  t*0  , 


(2.36) 


whereupon  (2.35)  becomes 


W^(eo(t))^0  ,  t*0  .  (2.37) 

This  In  turn  will  be  true  If 

W^(e)^0  ,  0se<61  .  (2.38) 

Materials  for  which  (2.38)  holds  are  said  to  be  non-hardening  in  the  first 
phase.  We  have  thus  shown  that  (2.33),  (2.36)  and  (2.38)  are  sufficient 
to  guarantee  that  c ( eQ ( x ) )  Is  a  monotonlcally  decreasing  function  of 
t,  so  that  (2.28)  Is  uniquely  Invertible  on  Rp. 

Finally,  the  rays  C  form  an  unambiguous  cover  for  all  x^O, 
tiO,  provided 

R={(x,t)|  xsO,  tiO],  RgORj.*!!)  (the  null  set).  (2.39) 

By  (2.31)  and  the  remarks  Immediately  following  it,  (2.39)  will  hold 
provided  c(eQ(0) )  * c(0) ,  which  In  turn  will  hold  If 

e  (0)  -  o  (0)  -  0  .  (2.40) 

o  o 

Physically,  (2.40)  Is  a  condition  that  the  Initial  loading  be  nonlmpul- 
slve.  In  summary,  conditions  sufficient  to  ensure  that  a  unique  ray 
C_  passes  through  each  (x,t)€R  when  0se<^  is  that  (2.33),  (2.36). 
(2.38)  and  (2.40)  are  satisfied.  We  shall  refer  to  these  conditions  as 
those  for  smooth  monotonic  loading  of  a  bar  whose  material  is  non-hardening 
In  the  first  phase. 
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D.  Discussion 

For  completeness  we  shall  briefly  discuss  the  phenomena  associated 
with  the  violation  of  any  one  of  conditions  (2.33),  (2.36),  (2.38)  or 
(2.40),  This  is  Instructive  not  only  for  Its  own  sake,  but  will  also  lead 
to  a  better  understanding  of  the  dynamic  phase  changes  associated  with 
the  violation  of  (2.27),  which  we  will  soon  consider. 

We  have  seen  that  a  unique  ray  C_  passes  through  any  point  (x,t) 
in  ftp  Simply  put,  the  Importance  of  (2.33),  (2.36)  and  (2.38)  lies 
solely  in  that  they  are  sufficient  to  ensure  that  c(eQ(x))  Is  a  con¬ 
tinuous  and  monotonically  decreasing  function  of  t,  which  by  (2.28) 
guarantees  that  a  unique  ray  C_  passes  through  each  point  (x,t)  in 
Rjj.  It  Is  then  (2.40)  which  provides  that  any  point  (x,t)  in  ft  is 
either  In  ftp  or  ftp  but  not  both. 

Condition  (2.33)  ensures  the  differentiability  of  c(eQ(x))  and 
so  reduces  the  investigation  of  the  monotonicity  condition  to  an  examina¬ 
tion  of  the  signs  of  derivatives.  Consequently  one  can  sacrifice  (2.33) 
provided  c(eQ(x))  remains  monotonically  decreasing.  Of  course  dis¬ 
carding  (2.33)  would  affect  the  differentiability  of  the  solution  (2.30). 

We  cannot,  however,  weaken  (2.33)  too  much  since  then  the  continuity 
of  c(e0(x))  is  threatened.  The  continuity  of  c(eQ(x))  can  be  at¬ 
tained  by  securing  the  continuity  of  the  constituent  functions  c  and  eQ. 
This  in  turn  is  guaranteed  by  requiring 


aQ(t)  €  C([0,co  ))  , 


W(e)  €CZ([0,61))  , 


(2.41) 
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the  second  of  which  was  an  original  requirement  on  the  energy  density 
(see  (2.14)).  Violations  of  either  of  (2.41)  are  not  necessarily  physi¬ 
cally  unreasonable;  discontinuities  In  o  (t)  are  associated  with  1m- 

o 

pulslve  loading,  while  Isolated  failures  of  (2.41 )^  Involve  piecewise 
differentiable  stress-strain  relations. 

Let  us  suppose  that  c(e0(x))  Is  monotonlcally  decreasing,  but 
that  (2.41)  falls,  so  that  c(cQ(x))  Is  discontinuous  for  some  t  =  td, 
say.  Then  the  characteristic  rays  abruptly  change  slope  at  x*0,  t  =  td, 
so  that  the  region  between  the  rays  parametrized  by  C^Ctj)  and  C_(td) 
consists  of  points  (x,t)  for  which  (2.28)  cannot  be  solved  for  x(x,t). 
Hence  we  must  extend  the  partial  solution  (2.30)  to  the  wedge-shaped 
region  enclosed  between  C_ ( t^)  and  Cft*).  Consider  first  the  case 
where  W*  (e)  has  an  Isolated  discontinuity  at  eaed-  Then  the  charac¬ 
teristic  rays  abruptly  change  slope  at  x * td  where  oQ(td)  ■W' (ed). 

On  both  C_(td)  and  C_(td),  (2.30)  Implies  that  e*eQ(td)  and 
v  *  <p( cQ( td ) ) .  Hence  the  solution  Is  completed  by  defining  e  =  eo(td) 
and  v  *  <p( eQ ( td ) )  In  the  region  enclosed  between  C_(tj)  and  C_ ( t^) . 

It  may  be  mentioned  that  this  solution  for  a  kinked  stress-strain  curve 
can  be  generated  as  the  limit  of  solutions  with  smooth  stress-strain 
curves  whose  slope  rapidly  changes  over  a  small  strain  interval,  Fig.  3. 

Now  we  consider  an  isolated  discontinuity  in  aQ(t)  at  t  =  td>  This 
prompts  a  discontinuity  In  eQ(t)  at  tatd>  Consequently  e  and  v  have 
different  constant  values  on  C_(td)  and  C_(td).  To  extend  (2.30)  into 
the  region  enclosed  between  these  rays  one  may  consider  end-tractions 
which  are  continuous  but  have  large  derivatives  on  some  finite  time 
Interval.  By  computing  the  solution  via  (2.30)  and  taking  the  limit  as 


the  Interval -length  vanishes.  It  may  be  verified  that  the  traction-dis¬ 
continuity  leads  to  the  well-known  centered  simple  waves  [16],  Fig.  4. 
Similarly  a  failure  of  (2.40)  also  gives  rise  to  a  centered  simple  wave. 

In  these  traction-impulse  constructions,  a  subtle  but  Important  point  Is 

that  If  W'(e)  Is  linear  between  e=e(t^)  and  e=e(t^)  -  i.e.,  in  (2.38), 
W'"(e)  =  0  for  eQ(tj)  ^eQ(tj)  -  then  all  the  characteristic  rays  in 
the  centered  simple  wave  have  the  same  slope.  Accordingly  the  simple  wave 
degenerates  to  the  single  characteristic  C_(td),  across  which  v  and  £ 
are  discontinuous.  Courant  and  Friedrichs  [16]  remark*  that  such  a  degen¬ 
erate  simple  wave  "...does  not  deserve  the  name  'shock'...".  Nevertheless, 

It  may  be  verified  that  the  shock  conditions  (2.12)  are  satisfied  across 
such  a  degenerate  simple  wave.  It  Is,  however,  a  particularly  well-behaved 

shock  In  that,  as  the  limit  of  classically  smooth  solutions.  It  Is  the 
limit  of  solutions  for  which  (2.23)  Is  satisfied.  Accordingly  (2.23) 
remains  satisfied  on  both  sides  of  such  a  shock,  so  preserving  straight 
line  characteristics  throughout  Rp,  each  of  which  still  propagates  a 
constant  value  of  v  and  e.  We  shall  see  that  such  a  fortuitous  situation 
Is  not  maintained  In  either  the  conventional  shocks  we  now  consider,  nor 
In  the  propagating  phase  boundary. 

We  have  so  far  examined  conditions  In  which  the  constructed  solution 
(2.30)  leads  to  regions  of  ft  devoid  of  any  rays  C_.  A  contrary 
state  of  affairs  occurs  when  c(eQ(t))  Is  increasing  on  some  Interval 
In  t.  Then  (2.28)  has  multiple  solutions  In  some  region  of  ftp.  This 
can  occur  if  either  c(e)  Is  increasing  on  some  Interval  or  if  aQ(t) 


Is  decreasing  on  some  Interval.  When  (2.33)  holds,  these  eventualities 
are  associated  with  a  violation  of  (2.38)  and  (2.36),  respectively.  In 
either  case,  the  result  Is  that  the  rays  C_  may  Intersect  each  other. 

Not  surprisingly.  In  those  regions  of  ftp  where  (2.28)  has  multiple 
solutions,  the  remedy  Is  to  discard  all  but  one  such  solution.  Presumably 
this  can  be  accomplished  by  Introducing  within  R  a  finite  number  of 
conventional  shock  curves  x*s^(t)  separating  regions  with  classically 
smooth  solutions  to  (2.8).  Across  these  curves  c  and  v  have  jump 
discontinuities  restricted  by  (2.12). 

We  note  that  In  the  x-t  plane,  the  first  conventional  shock 
x»s^(t)  so  generated  must  lie  to  the  right  of  the  leftmost  envelope  of 
intersecting  characteristics  C_.  The  solution  (2.30)  will  continue  to 
hold  for  x<s1(t)  provided  r(x,t)  is  taken  to  be  the  minimum  positive 
root  of  (2.28).  In  x>s1(t)  the  argument  leading  to  (2.23)  is  no  longer 
valid  since  we  cannot  Integrate  (2.21)  across  s^t)  without  using  (2.12) 
to  match  integration  constants.  Suppose  we  introduce  N  shock  curves. 
Then  all  we  can  conclude  from  (2.21)  Is  that  there  must  exist  2N 
functions,  K^,  K3***Kn+-|»  . . .Kjj+1 ,  such  ^ 

s1-l(t)  <  x  <  sj(t),  1*1,...,  N;  or  Sj(t)<x<0,  1*N,  (2.42a) 

then  C+  is  constant  on  curves  obeying 

^■-c(e),  v-<|.(e)  -k|(C+),  (2.42b) 

and  C  is  constant  on  curves  obeying 
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|r  -  -c(e) ,  v  +  <j>(e)  -  Kj(C_)  .  (2.42c) 

It  must  be  remembered  that  the  number  and  location  of  these  conventional 
shocks  x»s.j{t)  are  part  of  the  unknowns  in  the  problem.  After  the 
introduction  of  these  shock  curves,  any  point  (x,t)  €R  should  lie  on 
the  intersection  of  unique  C_  and  C+  characteristic  curves.  Moreover, 
as  the  shocks  were  introduced  to  eliminate  multiple  characteristics  C_ 
through  certain  points  In  R,  we  expect  that  the  shocks  themselves 
should  not  give  rise  to  additional  characteristics  C_.  In  other  words 
all  curves  C_  should  intersect  either  t«0  or  x»0.  Indeed  the 
requirement  that  all  curves  C_.  should  originate  on  initial  or  boundary 
data  is  a  special  case  of  an  admissibil ity-or  entropy -condition  formu¬ 
lated  by  Lax  for  more  general  hyperbolic  systems  [18], 

So  far  we  have  dealt  exclusively  with  the  soft  device  governed  by 
the  traction  boundary  condition  (2.11).  The  results  for  the  hard  device 
governed  by  the  kinematic  boundary  condition  (2.10)  are  also  quite  similar 


for  the  case  where  e(x,t)<B1  .  For  convenience  we  define 

vQ(t)  »  0Q(t)  ,  t*t0  .  (2.43) 

Then  smooth  solutions  to  (2.8),  (2.9)  and  (2.10)  are  given  by 

etO,  v  =  0  for  x  <  -c($(vQ(0)))t  ,  (2.44) 

which  is  the  tranquil  region  in  front  of  the  leading  disturbance.  While 
v=vo(x(x,t))  ,  e  =<l>(vo(T(x,t)))  ,  (2.45) 


in  the  disturbed  region 


xi-c(f(vo(0)))t  .  (2.46) 

Here  r(x,t)  is  given  implicitly  by 

x * -c($(vQ(T)))(t  -  t)  .  (2.47) 

The  stipulation  that  e(x,t)  is  less  than  the  value  B-j  will  hold  provided 

vQ(t)  <  d>(61 }  ,  t*0  ,  (2.48) 

while  (2.47)  will  have  a  unique  solution  in  the  disturbed  region  provided 
both 


vQ(0)«0  ,  (2.49) 

and  c($(vo(t)))  is  a  continuous  and  monotonically  decreasing  function 
of  t  .  If  we  consider  materials  which  are  nonhardening  in  the  first 
phase —so  that  (2.38)  holds— then  a  condition  sufficient  to  ensure  the 
above  monotonicity  condition  is  that  vQ(t)  be  differentiable  and 

v  (t)*0  .  (2.50) 

o 

When  (2.38),  (2.49)  or  (2.50)  are  violated,  one  encounters  difficulties 
similar  to  those  of  the  specified-load  problem.  In  these  instances 
analogous  remedies  are  invoked. 

Finally,  we  remark  that  the  consideration  of  compressive  end-loadings 


from  an  initially  undisturbed  rest  state  would  require  the  stress  response 
to  be  known  for  negative  strains.  If 
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W*(e)>0  ,  |T(e)sO  for  e<0,  (2.51) 

the  associated  dynamical  fields  are  smooth  provided  there  Is  no  unloading. 
In  contrast  to  the  behavior  encountered  In  the  anal agous  piston  problem 
of  gas  dynamics  [16]. 


3.  Characteristic  Theory  when  Two  Phases  are  Present. 

So  far  we  have  only  considered  loadings  which  permit  the  material 
to  remain  In  the  first  phase  throughout  the  bar  for  all  time.  We  now 
consider  the  consequences  of  end  conditions  that  lead  to  e  >  B-j  *  so 
that  solutions  are  not  restricted  to  the  Initial  hyperbolic  regime  of 
strain 

A.  Formulation  of  the  Two-Phase  Problem  for  the  Soft  Loading  Device. 

We  now  return  to  the  traction  problem  governed  by  (2.8),  (2.9), 
(2.11)  and  (2.12).  As  In  the  situations  leading  to  conventional  shocks, 
phase  boundaries  arise  even  with  arbitrarily  smooth  data.  To  focus 
attention  on  the  specific  phenomenon  of  phase  changes,  we  shall  assume 

H(e)6rfD([0,m)),  <Jo(t)€(f“([0,<B)).  (3.1) 

As  In  the  previous  section  we  further  assume  that  (2.36),  (2.38)  and 
(2.40)  hold.  For  convenience  we  collect  these  together  below: 

1)  aQ(0)  *  0*  aQ(t)  tiO  , 

(3.2) 

11)  kT'teJsO  ,  0se<8i 

We  have  already  shown  that  If,  In  addition,  (2.27)  holds,  then  there 
Is  a  classically  smooth  solution  (2.30)  with  e  <  B-|  throughout  the  bar 
for  all  time.  Thus  we  only  need  consider  the  case  where  aQ(t)  attains 
or  exceeds  the  value  W*(B-|)  s  .  Let  tg  be  the  time  for  which 
(2.27)  Is  first  violated  so  that  t * tg  Is  the  minimum  value  of  t 
satisfying 


Now  (2.15)  still  guarantees  the  existence  of  a  unique  function  fcQ(t) 
defined  on  the  Interval  0<t<tg  such  that 

W' (eQ(t))  *  aQ(t)  ,  C 


OseQ(t)  <61  .  (3.5 

If  we  examine  the  previous  solution  (2.30).  with  eQ(t)  defined  through 
(3.4)  and  (3.5),  we  find  from  (2.18),  (2.15)  and  (3.3)  that  c(eQ(T(x,t) 
approaches  zero  as  t(x,t)  approaches  the  value  t0  .  As  a  result, 

P 

the  rays  C_  In  the  x-t  plane  become  Infinitely  steep  as  their 
t  -  Intercept  approaches  tg  .  Moreover  the  set  of  rays  C_  form  a 
nonintersecting  cover  of  (0,t)  |t  *tg};  see  Fig.  5.  Thus  if  we 

append  the  restriction  r<tg  to  (2.30),  we  obtain  a  smooth  solution 
everywhere  except  for  points  (x,t)  with  x*0  and  titg. 

It  Is  natural  to  ask  whether  the  set  of  such  points  with  x*0  and 
t^tg  could  form  some  sort  of  stationary  "end-shock".  For  this  to  be 
the  case  s(t)»0,  so  that  (2.12)g  and  (2.30)  would  now  imply 

W'UfO.t))  -  ¥*(eo(x(0",t)))  »  0-  ,  titg  .  (3.6 


This,  with  (2.11),  necessitates  that  (3.3)  must  hold  for  all  times 
t  greater  than  the  value  tfl.  One  also  sees  that  If  a  (t) 

Is  maintained  at  the  value  Og  for  some  Interval,  say  t^stst^,  then 
we  may  extend  (2.30)  smoothly  onto  x«0  and  tg^t^t^  by  e(0,t)  s  6 
and  v(0,t)  * ^(B-, )  for  tgSt*tK  •  We  are»  however,  faced  with  the 

problem  of  not  being  able  to  construct  a  solution  by  this  method  for  the 
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set  of  points  x«0  and  t  >  t^.  Consequently  such  stationary  end-shocks 
are  precluded  if  a  (t)  exceeds  the  value  aa.  For  simplicity  we  shall 
suppose  that  t*tg  is  the  only  root  of  (3.3).  Then  it  follows  from 
(2.11),  (2.15)  and  (3.2)1  that 

e(0,t)  >  B2  ,  t  >  tg  .  (3.7) 

Now  (3.4)  and  (3.5)  ensure  that  e(O,t)<0i  is  consistent  with  the 
loading  condition  (2.11)  for  t<t„.  The  relations  (2.15),  (2. 17)  and 

P 

(3.3)  show  that  e(0,tg)  may  take  on  either  the  value  B-j  or  the  value 
02*  Hence  it  is  natural  to  seek  solutions  where  e(0,t)  jumps  from 
6l  to  B2  at  0nce  one  concludes  that  the  induced  end-strain 

must  be  discontinuous,  it  Is  clear  that  a  jump  from  the  first  to  the 
second  ascending  branches  of  the  stress-strain  curve  could  occur  at 
the  end  of  the  bar  for  any  t  In  the  interval  where  (see 

Fig.  1), 


W  *  woo(V  •  (3-8) 

In  the  face  of  this  nonuniqueness,  we  opt  to  seek  solutions  for  which 
the  end-strain  lies  below  the  first  ascending  branch  for  as  long  as 
feasible,  as  there  is  no  mathematical  need  to  introduce  a  second  phase 
until  the  loss  of  hyperbolicity  is  Imminent  at  t  =  tg  .  Consequently  we 
shall  define  a  particular  inverse  for  W/(e),  denoted  by  T(»).  such 


that 


(3.9) 


T:[0,cd  )  [0,6-j ]  U  (02»®  )  » 

T(W'U)}  *  e  for  >  £>62* 

By  (3.1),  T  Is  differentiable  to  all  orders  except  at  where  T  is 
discontinuous  (Fig.  6).  With  this  choice  of  inverse,  the  induced  end- 
strain  is  given  by 


eQ(t)  5  T(ao(t))  ,  t*0  .  (3.10) 

We  note  that  e  (t)  Is  discontinuous  at  t*tfl.‘  For  all  other  t^O, 
eQ(t)  is  differentiable  to  all  orders,  with  the  first  derivative  again 
given  by  (2.32).  Moreover,  since  W*(6, )«0  and  W"(62)>0  It  follows 
that  the  right  hand  limit  of  all  derivatives  of  eQ(t)  exist  at  t«tg, 
whereas  these  derivatives  become  unbounded  as  t  approaches  tfl  from  the  left. 

P 

Since  a  stationary  end-shock  Is  Incompatible  with  an  end-load  greater 
than  the  value  a^.  It  follows  that  the  high-strain  phase  will  not  be 
confined  to  the  end  of  the  bar,  but  that  a  phase  boundary  will  emerge  at 
the  end  of  thebarattlme  t_  and  subsequently  travel  Into  the  Interior. 

P 

In  front  of  this  (unknown)  phase  boundary,  which  wedenoteby  x*s(t), 
the  strain  and  velocity  fields  are  still  given  by  (2.30),  while  for 
s(t)<xs;0  we  seek  solutions  with  e  >  62  • 

Note  that  for  our  particular  choice  of  Inverse  to  the  stress  response, 
the  phase  boundary  Is  not_  genearated  at_  the  Maxwell  stress  a^,  which  In 
the  equilibrium  case  Is  the  stress  where  the  configuration  loses  absolute 
stability.  Instead  the  phase  boundary  emerges  when  the  end-load  reaches 
a„,  which  in  the  equilibrium  setting  is  the  stress  at  which  the  config- 
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uratlon  loses  metastability.  Consequently  we  are  assuming  a  type 
of  "superstraining"  analogous  to  the  superheating  a  liquid  undergoes 
when  heated  above  Its  boiling  point.  This  would  seem  reasonable  for  a 
smoothly  applied  end-loading;  In  this  case  we  regard  the  phase-boundary 
as  emerging  in  response  to  the  Incipient  loss  of  hyperboll city. 

We  now  extend  the  definitions  (2.18)  and  (2.19)  to  the  high-strain 
phase  by  writing 


Once  again  <j>(e)  Is  monotonic  on  Its  domain  of  definition,  so  that  we  can  ex¬ 
tend  the  previously  defined  Inverse  9  to  the  range  of  <J>  as  given  by 
(3.12).  <f(z)  Is  continuous  on  Its  domain  0s;z<11m  <J>(e),  and  by  (3.1) 

G-*^D 

is  differentiable  to  all  orders  except  at  z  =  4>(6i )  (=  <p(a2) ) •  At 

za<J>(Bi)»  the  one-sided  derivatives  of  4>(z)  exist  from  the  right  but 
are  unbounded  from  the  left,  ^'(z)  and  the  higher  derivatives  may  still 
be  found  from  (2.20) 

Denote  by  R|  and  the  regions  in  the  low-  and  high-strain  phases 
respectively,  so  ttiat 
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ftj  ■  {(x,t)|(x*0,  titg)  U  (x  <  s(t),  t>tg)3  , 
<*2  ■  £(x .t) | s( t )  <x*0  ,  t>tg3  • 


1 


(3.13) 


We  shall  refer  to  fej  as  the  region  in  front  of  the  phase  boundary  and 
ft2  as  the  region  behind  the  phase  boundary. 

In  ftp  the  dynamical  fields  are  still  given  by  (2.30)  with  the 
restriction  that  T<tg.  In  crossing  the  phase  boundary,  we  require 
that  the  shock  conditions  (2.12)  be  satisfied.  In  ftg,  the  extension 
of  c  and  4  Into  the  second  hyperbolic  regime  of  strain  allows  us  to 
adopt  a  scheme  like  that  discussed  In  the  preceding  section  for  finding 
solutions  behind  conventional  shocks. 

In  we  sha11  parametrize  the  positively-sloped  characteristic 
curves  by  the  time  at  which  they  Intersect  x»0.  This  Is  accomplished 
by  letting  T+(x,t)  be  the  time  at  which  the  positively-sloped  character¬ 
istic  curve  passing  through  the  point  (x,t)  In  ^  terminates  on  x»0. 
Notice  that  from  this  definition  follows  T+(0,t)  -  t.  Similarly 
let  t"(x,t)  parametrize  the  negatively-sloped  characteristic. curve  In 
which  originates  on  x-0  at  t»x.  Under  the  extended  definition 
of  c(e),  (2.21)  continues  to  hold  In  R^.  As. before,  these  equations 
may  be  Integrated  to  yield  the  Rlemann  Invariants  <2  and  ,  where 

v(x,t)  -  <p(e(x ,t) )  ■  K2(x+(x,t)),  (3.14a) 

on  curves  parametrized  by  T+(x,t),  for  which 


■  c(e(x,t))  ,  T+(0,t)  »  t  . 


dx 

3F 


(3.14b) 


Similarly, 

v(x,t)  +4>(e(x,t))  *  ^(r’Cx.t))  ,  (3.14c) 

on  curves  parametrized  by  -r’(x,t),  for  which 

-  -c(e(x,t))  ,  T_(0,t)  3  t  .  (3.1 4d ) 

Here 

T+:R2^[tg,a)  )  ,  T'cRg-^Ctg.oo  )  ,  (3.1 4e ) 

K^Ctg.oo^R  ,  K’:[te,ffl  )  +  R  ,  (3.14f) 

constitute  unknowns  to  be  determined.  There  Is  a  final  unknown,  namely 
the  phase  boundary  location  s(t).  We  have  seen  that 

s(t)iO  for  titg  ,  s(tg)  =  0  •  (3.15) 

We  shall  seek  s(t)  GC*  ([tg,a>  ))  flC®  ((tg.oo  ))  .  The  shock  conditions 
(2.12)  become 

[€(s(t),t)]5(t)*[v(5(t),t)]-0,  [v(s(t),t)]5(t)+I[U'(e(s(t),t))]-0 

for  tzta  .  (3.16) 

Here  e(s(t)",t)  and  v(s(t)’,t)  are  given  by  (2.30)  with  x(s(t) ,t)<tg. 
We  must  also  satisfy  the  boundary  condition 

e(0,t)  » eQ(t)  ,  t*tg  .  (3.17) 

For  (x,t)  In  r2,  (3.14)  shows  that  the  velocity  and  strain  are  given 


(3.18) 


v(x,t)  *  jKg  (x  (x,t)) +j-K2  (t  (x,t)), 
e(x,t)  -^Q-Kg  (x”(x,t))  -  j-Kg  (x+(x,t))j  . 

We  may  use  (3.14)  and  (3.17)  to  eliminate  as  follows: 

2<fr(e0(t))  -  2<|>(e(0,t))  *  K"(T_(0,t))  -  K2(x+(0,t))  -  K^(t)  -  K^t),  (3.19) 

so  that 

K“(t)  ■  Kg(t)  +  2 <t> ( eQ ( t ) )  .  '  (3.20) 

Thus  (3.18)  becomes 

v(x,t)  ■  jK*  (x+(x,t))  +  jK*  (x"(x,t))  +4>(E0(x"(x,t))):, 

(3.21) 

e(x,t)  K2(x+(x,t))  +  j  K2(x"(x,t))  +4»(eo(x"(x,t)))j  . 

Collecting  our  results  so  far,  we  may  summarize  the  problem  as  follows: 
When  W(e)  and  aQ(t)  are  subject  to  (2.14),  (2.15),  (2.16),  (2.17), 
(3.1),  (3.2),  and  with  EQ(t),  c(e),  $(e),  $(z)  defined  through  (3.10), 

(3.11),  (3.12),  we  seek  four  functions  K2(z),  x+(x,t),  x"(x,t)  and  s(t), 
such  that: 

(i)  s(t)  ec1  ([tg.co ))  nc®  ((tg,oo )),  s(t)  so,  s(t0)*o 

and  s(t)  partitions  R  Into  Rj  and  R2  as  defined  by  (3.13); 

(11)  In  r.j  ,  e(x,t)  and  v(x,t)  are  given  by  (2.30)  with  the  restric¬ 
tion  x<t  ; 


(111)  In  r?,  e(x,t)  and  v(x,t)  are  given  by  (3.21); 


( 1  v)  T+:  Rg^Ctg.oo )  ,  T+(0,t)*t  and  t+ 
for  which  ■  c(e(x,t)); 

(v)  T’l.  ftg-^Ctg.OD  )  ,  x“(0,t)  a  t  and  t" 

for  which  ^  ■  -c(e(x,t)); 


Is  constant  on  curves 


Is  constant  on  curves 


(vl)  The  jump  conditions  (3.16)  are  satisfied. 


(3.22) 


Solutions  to  (3.22)  give  a  single  phase  boundary  preceding  a  region 
of  high-strain  material  phase  and  propagating  into  a  region  of  low-strain 
material  phase.  In  both  of  these  regions  the  strain  and  velocity  fields 
are  classically  smooth  solutions  to  (2.8).  It  should  be  emphasized 
that  the  existence  of  such  a  smooth  solution  behind  the  phase  boundary  is 
by  no  means  guaranteed.  In  this  regard  recall  that  the  material  condition 
V"'  (e)  <0  for  0<e<8|,  Is  necessary  for  the  existence  of  smooth  fields  in 
front  of  the  phase  boundary,  whereas  we  have  not  as  yet  Imposed  a  comparable 
condlton  on  W"'(e)  for  e^o^. 

B.  Statement  of  the  Two-Phase  Problem  for  the  Hard  Loading  Device. 

Before  investigating  (3.22)  further  we  shall  give  a  statement  of 
the  corresponding  free  boundary  problem  for  the  hard  device  governed  by 
(2.10).  Here  we  again  define  c(e),  $(e)  and  $>(z)  by  (3.11)  and  (3.12). 
we  consider  a  non-lmpulslve,  smooth  loading  of  a  bar  of  material  that  is 
non-hardening  in  the  first  phase,  we  are  assured  that  (2.38),  (2.49)  and 
(2.50)  hold.  It  Is  then  possible  for  the  bar  to  remain  everywhere  in 
the  first  phase  until  time  tfl  ,  where  t * tfi  Is  the  largest  root  of 


vrt(t)  =  <f(Bn)  . 


(3.2 
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We  shall  only  consider  the  case  In  which  (3.23)  has  a  unique  root,  rather 
than  an  Interval  of  roots.  For  t  >  tg  near  the  end  x»0  It  Is  neces¬ 
sary  that  the  material  be  In  the  high-strain  phase.  If  x  =  s(t) 
denotes  the  resulting  propagating  phase  boundary,  we  seek  three  addition¬ 
al  functions  K^z),  T+(x,t),  and  T’(x,t)  such  that: 

(1)  s(t)  ec^Ctg.a)  ))  nc°°((tg,a>  )),  s(t)  sO,  s(tg)*0  and  s(t) 
partitions  ft  Into  ft^  and  as  defined  by  (3.13); 

(11)  In  ,  e(x,t)  and  v(x,t)  are  given  by  (2.44),  (2.45),  (2.46)  and 

(2.47)  with  the  restriction  r<tg  ; 

(111)  In  ft2,  e(x,t)  and  v(x,t)  are  given  by 

v(x,t)  -  (t+(x,£))  -  (t’(x,t))  +vo(x"(x,t)), 

e(x,t)  *$(-7*2  (T+(x.t))  *7K2  (T"(x,t))  +v0^t^ ’(x,t))); 


tg.ao  ), 

T  +  (0,t)  *  t 

and  x+ 

Is  constant  on  curves 

for  which  ^  ■ 

c(c(x,t)) ; 

T’rR^Ctg,®  ), 

T_(0,t)  *t 

and  t” 

Is  constant  on  curves 

for  which  ^  * 

-c(e(x,t)); 

(vl)  The  jump  conditions  (3.16)  are  satisfied. 


(3.24) 
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4.  A  Special  Class  of  Materials. 

The  application  of  characteristic  theory  leading  to  the  formulation 
(3.22)  requires  that  a  unique  member  of  both  families  of  characteristic 
curves  should  pass  through  each  point  (x,t)  In  R.  For  portions  of  the 
bar  In  the  low- strain  phase,  we  were  able  to  use  the  Initial  conditions 
(2.9)  to  show  that  the  Rlemann  Invariant  associated  with  the  positively- 
sloped  family  of  characteristic  curves  vanishes  Identically  (see  the 
discussion  preceding  (2.23)).  However.lt  Is  necessary  to  stipulate  that 
W*"(e)sO  for  Os  e<S^  In  order  to  guarantee  that  the  negatively-sloped 
family  of  characteristic  curves  (here  rays) forms  an  unambiguous  cover  of 
fcp  We  now  Inquire  as  to  what  requirement  (If  any)  on  W(e)  and  Its  derivatives 
In  Ocl^  will  guarantee  that  each  of  the  two  families  of  characteris¬ 
tic  curves  will  unambiguously  cover  R2*  Should  these  covering  proper¬ 
ties  fall,  x+  or  t*  would  either  be  undefined  or  multi-valued  at 
certain  locations  In  R2,  and  the  problem  governed,  by  (3.22)  would  not 
be  well -posed.  For  example,  multi-valuedness  of  t+  or  x”  would  sug¬ 
gest  the  presence  of  conventional  shocks  In  the  high-strain  phase  portion 
of  the  bar.  Unfortunately  we  do  not  know  of  a-priori  conditions  that 
are  both  necessary  and  sufficient  for  x+(x,t)  and  x"(x,t)  to  be 
both  defined  and  single-valued  for  all  (x.tJSRg.  However, in  what  follows 
we  give  a  condition  on  the  material  sufficient  to  ensure  the  existence 
and  single-valuedness  of  x+  and  x”  throughtout  Rg. 

We  shall  henceforth  consider  only  materials  with  the  property  that 
there  exists  a  strain  <5,  a2<<$<62,  such  that 


W4"  (e)  =  0  for  eifi  . 


(4.1) 
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Thus  there  are  constants  E  >  0  and  D  such  that 

M#(e)*Ee  +  D  for  e^6  .  (4.2) 

In  other  words,  we  shall  deal  only  with  energy  densities  for  which  the 
stress  Is  a  linear  function  of  strain  for  values  of  strain  greater  than 
6,  where  6  lies  below  the  second  ascending  branch,  but  is  less  than 
82-  A  diagram  of  such  a  stress-strain  curve  satisfying  (2.14),  (2.15), 
(2.16),  (2.17),  (2.38),  (4.1)  and  (4.2)  Is  found  In  Fig.  7. 

A.  Reduction  of  the  Problem  to  a  Pair  of  Differential -Del ay  Equations. 
From  (4.2),  (3.11)  and(3.12)  It  then  follows  that  for  such  a  material 


c(e)  -  yj ,  £*6  , 

(4.3) 

<J >(e)  *  0(6) +^^"(£-6),  e*6  , 

(4.4) 

$(z)  »  <J>(6)  +  5»  z  ^<p( 6)  . 

(4.5) 

Turning  our  attention  to  the  hard-device  problem  in  the  form  (3.22),  we 
have  by  (4.2),  (3.9)  and  (3.10)  that 

eQ(t)  *  ^(t)  -D)  >B2  for  t  >  tg  .  (4.6) 

We  seek  solutions  which  are  In  the  linear  range  corresponding  to  the 
high-strain  phase,  i.e.  c>6.  The  constant  acoustic  speed  (4.3)  in 
ft2  then  allow  us  to  integrate  ( 3 . 22) ^ y  and  (3.22)y  and  solve  for 
x+(x,t)  and  x”(x,t).  The  result  Is 


t  (x,t)  *  t  -  x  ,  j 

•  nr  r 

T  (x,t)  »  t +  £  X 

\ 

J 

Thus  for  (x.tJeRg  ,  (3.21),  (4.4),  (4.5)  and  (4.7)  yield 
*(x,t)  4k*  (t-^f  x  )  4k*  (f  Jfx)  + 

c(x.t)  -  -J-  Jfx)*  €0(t+^fx). 

It  is  convenient  to  define  the  new  function  f  in  terms  of  K*  by 

f(z )  * 7 y/IK2  ^Fz )  *  76  +  7  5) » 

whereupon  (4.8)  becomes 

v<x.t).  jFfu-vi^vifi-x-VTt)*^^*). 

e(x,t)  - -f(x-JJt)  *  f(-x-JJt)  ♦  e0(t+JJx)  ■ 


That  e(x,t)  and  v(x,t)  are  each  linear  combinations  of  functions 
of  the  variables  x  +  and  x-  yj^ t  is  a  consequence  of  the  fact 
that,  for  the  linear  portion  of  the  stress-strain  relation,  (2.8)  becomes 


2  2 
3  u  E  3  u  _ 


3?  p  3? 


=  0  . 


Differentiating  (4.11)  one  shows  that  s  and  v  aiso  satisfy  (4.T1).  Con¬ 
sequently  we  have  simply  retrieved  a  form  of  D'Alemberts  solution  to  the 

i 

classical  wave  equation.  Note  also  that  in  ft2  both  families  of 


characteristic  curves  are  straight  rays  with  the  respective  slopes 

.  Here,  however,  neither  family  of  characteristic  rays 
propagate  constant  values  of  e  or  v. 

Before  we  describe  how  the  results  of  this  analysis  simplify  the 
formulation  (3.22),  we  make  a  few  observations.  From  the  jump  conditions 
(3.16)  we  see  that  the  square  of  the  phase  boundary  velocity  is  given  by 


s2(t) 


1  [H* (e(s(t) ,t))1 
P  [e(s(t),t)] 


(4.12) 


and  so  for  a  normalized  density  (p*l)  is  given  by  the  slope  of  the 
secant  line  to  the  stress-strain  curve  between  e(s(t)+,t)  and 
e(s(t)~,t).  On  the  other  hand,  by  (3.11),  the  square  of  the  acoustic 
speed  at  any  location  is 

c2(e(x,t))  -  iu*  U(x,t)),  (4.13) 

2 

and  so, for  a  normalized  density,  c  (e)  is  given  by  the  slope  of  the 
stress-strain  curve  at  e(x,.t). 

In  the  problem  we  are  considering,  the  stress-strain  curve  for  a 
strain  associated  with  ft2  has  slope  E,  which  is  necessarily  greater 
than  the  slope  of  a  secant  line  projected  back  to  a  point  on  the  first 
branch  when  e(s(t)+,t)  is  near  ^  and  e(s(t)~,t)  is  near  8-j .  Hence 
(4.12)  and  (4.13)  Imply  that  for  some  A>0, 

| s ( t) I  <  c(c(s(t)+,t)  tg^tstg  +  A  •  (4.14) 


For  the  above  time  interval  the  phase  boundary  Is  subsonic  with  respect  to 
the  material  behind  it.  In  this  investigation  we  shall  restrict  attention 


to  phase  boundaries  for  which  this  is  true  for  all  time.  Hence  we  impose 


the  condition 


|s(t)|  <  c(e(s(t)+,t)  * 


'6  .  » 


(4.i5: 


an  Important  consequence  of  this  restriction  being  the  lower  bound 


t*tji . 


(4.16! 


Integrating  (4.16)  and  using  the  conditions  (3.15),  we  arrive  at 


v 0 


(4.17 


By  (3.13),  we  further  have 


for  (x,t)  €R2 


Hence  f(z)  need  only  be  defined  on  z* 


(4.18 


*  yJTh  * 


which  we  shall  hence¬ 


forth  take  to  be  the  domain  of  defl niton  for  f.  This  could  have  been 
predicted  from  (4.9)  and  (3.1 4f ) .  Notice  also  that  the  bound  (4.16) 
assures  that,  for  Increasing  t,  s(t)  Intersects  members  of  the  negatively 
sloped  family  of  characteristic  rays  In  ft2  which  arise  progressively 
later  on  x  a  0,  tit-. 

P 

We  can  Incorporate  the  Initial  conditions  (2.9)  into  the  resulting 
end-strain  by  setting  eQ(t)  equal  to  zero  for  t<0.  Then  by  (3.10) 
and  (4.6)  _ 


T(o0(t)) 


0  i  t  ^  tg , 


l{oJt)-D)  t>t- 


(4.19 


Note  that  e  (t)  Is  continuous  everywhere  except  at  t  ■  te  and  con- 

tinuously  differentiable  to  all  orders  except  at  t  *  0  and  t  =  t^.  We 

also  define  a  (t)*0  for  t<0  so  that  (3.10)  holds  for  all  time, 
o 

Finally  we  define  x(x,t)“0  for  (x,t)€Ry  .  7  These  extensions  allow 
us  to  use  the  same  formulation  In  both  the  disturbed  and  undisturbed 
regions  of  the  bar. 

By  means  of  (4.10)  the  problem  formulated  In  (3.22)  Is  reduced  to 
finding  functions  s(t)  and  f(z)  for  which 


(1)  s(t)  CC^Ctg,®  ))  nc°°  ((tg,co  )),  s(t)-i 0  and 


s(tg,)  *0; 

(4. 

(11)  f(z):  (-a),  -  ^|t6]-*-F  ; 

(4. 

(111)  [e(s(t),t)]s(t)+[v(s(t),t)]-0  for  t>tg. 

(4. 

and 

Cv(s(t),t)]s(t)*I[W#(e(s(t).t))]-0  for  t>t&. 

(4. 

where 

[v]  =  Cv(s(t),t)] 


•Jff  >/f *yff- >/f  *>  *>/F £o(t  *(eo(T))> 


(4.: 


[e]  =  [e(s(t),t)]  * 

-f(s(t) /It,  +  f(-s(t)-.  /ft,  +  e  (tt  If 5(t»-€(t),  (4.: 


[tf]  =  W(c(s(t),t))]  * 


-Ef(s(t)  -  +  Ef(-s{t)  -yfft)  +  Ee0(t  +Jfs(t))  +  D-M»  U0(t))  ,  (4.20g) 

and  x*t(s(t),t)  Is  the  unique  root  less  than  tg  of 

s(t)  -  -c(e0(T))(t-Tj  .  (4.20h) 

In  summary,  for  an  Induced  end-strain  given  by  (3.10)  the  problem 
formulated  In  (4.20)  completely  characterizes  the  phase  boundary 
and  dynamical  fields  associated  with  the  smooth  monotonic  loading 
(In  a  soft  device)  of  a  bar  of  material  that  Is  non-hardening  In  the 
first  phase  and  satisfies  (4.2).  The  system  (4.20)  consists  of  two 
differential -delay  equations  for  the  unknown  functions  f(z)  and  s(t). 
Notice  that  s(t)  appears  In  the  retarded  argument  of  f  so  that  the 
amount  of  delay  constitutes  part  of  the  solution. 

It  will  be  convenient  to  consider  the  following  linear  combinations 
of  (4.20c)  and  (4.20d): 

+  »)[*]-  •  <4-21b 
Upon  substitution  from  (4. 20e,f,g)  Into  (4.21),  we  arrive  at 

2^I(s(t)  +>/|)f(-s(t)  -^t)  +  r1(s(t),s(t),t)  =  0  ,  (4.22a) 


and 


-  •  •  ,»V,  «  ,  *  * 

*  *  4  * -r  *  *  , 
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+^f)f(s(t)  -^|t)  +  r2(l(t),s(t),t)  ■  o 


t  ^  t- 


(4.22b) 


where 


)<j»(eo(T(s,t)))-la0(x(s,t))  , 


(4.23a) 


and 


r2(s,s,t)  -^|"se0(T(s,t))  +  (s  -^|)^(e0(T(s,t))) 


+  £  a0(x(s,t)). 


(4.23b) 


We  shall  use  (4.22)  and  (4.23)  In  the  forthcoming  asymptotic  analysis. 

B.  An  Alternative  Formulation  Involving  an  Integral  Equation. 

It  is  also  possible  through  (4.22)  and  (4.23)  to  eliminate  f  by 
recasting  the  problem  as  an  Integral  equation  for  s(t). 

Define 

t 

Y(t)  -  j  f(-J|e)da  for  t  as  t0  .  (4.24) 

p 

Then  (4.22a)  can  be  written  as 


2i .^Y(^|  s(t) +t) +ri(s(t),s(t),t)  *  0  for  tatg  .  (4.25) 


Integrating  (4.25),  using  Y(tg)«0  and  s(tg)*0,  we  find 


(4.26) 


t 

|y(jfs(t)+t)  +  |  iysoo, 

trt 


s(a),a)daaO  for  tit 


6  * 


Similarly  (4.22b)  leads  to 


t 

2  py  (  _-jFs(t)  + 1)  +  |  r2(s(«),s(«),a)da  ■  0  for  tit^  .  (4.27) 

By  (4.15),  both  -^”s(t)  +t  and  ^jj^s(t)+t  are  continuous  and 

monotonlcally  Increasing.  Thus  for  any  titfl  ,  there  exist  times  t 

p  a 


and  tb  such  that 


a  ^(tb  )+tb  .  (4.28) 


Note  that 


taititbitB  ’ 


(4.29) 


Evaluating  (4.26)  at  t  and  (4.27)  at  t.  we  may  eliminate  Y(t)  to 

a  D 

obtain 

ta  tb 

f  r.(s(a),s(a),a)da- f  r?(s(a) ,s(a) ,a)da  = 0  .  (4.30) 

tg  tB 

To  simplify  (4.30)  we  note  that  by  (4.28)  and  (4.23a),  the  first  term  In 
the  Integral  on  the  left  Is 

ta  r  r-  1 

J  (s(a)  +y|-)o0(a  +^|-s(a))da  «  ^  J  aQ(a)da  .  (4.31) 

t  „  8 


With  the  help  of  (4.23),  (4.29)  and  (4.31),  (4.30)  yields  the  following 
equation  for  s(t): 


+  — —  o  (t(s ,t) )  ,  (4.3; 

2/pT  0 

and  t  ■  t  (t)  and  t.  *t.(t)  are  the  roots  of 
a  a  b  d 

•  (4-3 

Notice  that  we  can  evaluate  the  left-hand  side  of  (4.32a)  as  well  as  the 
first  terms  In  the  integrals  on  the  right-hand  side.  For  an  induced 
end-strain  given  by  (3.10),  system  (4.32)  also  completely  characterizes 
the  phase  boundary  and  dynamical  fields  associated  with  the  smooth 
monotonic  loading  (in  a  soft  device)  of  a  bar  of  material  that  is  non¬ 
hardening  in  the  first  phase  and  satisfies  (4.2).  Of  the  two 
formulations  (4.20)  and  (4.32),  we  shall  consider  only  (4.20)  in  the 
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5.  Emergence  of  the  Phase  Boundary. 

We  now  study  the  behavior  of  solutions  of  the  system  (4.20)  near 
t»tQ,  at  which  time  the  phase  boundary  first  appears.  Since 

P 

s(tg)  *  0,  we  have 


;(t8 tB  * 


(5.1) 


so  that  at  the  Initial  time  tg  ,  the  two  arguments  of  f  Involved  In 
(4.20)  are  the  same.  Differential -del ay  equations  whose  arguments 
coalesce  Initially  are  called  singular,  and  In  general  they  permit  the 
specification  of  only  the  Initial  values  of  the  unknown  functions;  see 
[22].  In  our  case,  t.he  Initial  value  of  s  Is  given:  s(tg)*0  .  More¬ 
over,  we  may  let  t  +  t*  In  (4.20)  to  determine  f (  -  ^|  t g).  Recalling 
that  eQ(t)  Is  discontinuous  at  t«tg  ,  we  have 

11m+  eQ(t>  /|s(t)) -eQ(tg) -B2  ,  (5.2) 

t  ^  tfl  * 


11m+  e0(T(s(t),t))  » eQ(tg)  *  81  ,  (5.3) 

t  +  tg 

so  that  (4.20e)-(4.20g)  »  when  evaluated  at  t»tg,  become 

[vV2^f(>/ite)VI62-*(6i)  •  (5-4) 

[e]  -B?  -6,  ,  (5.5) 

«  0  . 

*8 


(5.6) 


We  may  then  Infer  from  (4.20c)  and  ( 4 . 2 Od )  that 

<B2  -  B1  >5UB)  *  *(8, )  -  0  ,  (5.7 

(2^5ft*^IV*>/FB2-*(6l))5(tB)'0  '  <5'8 

These  equations  Imply  that 

s(tg)  ■  0  ,  (5.9 

f< -VFv  *?(>/!  *(ei>- O'  (5-,° 

We  observe  that,  from  (5.4)  and  (5.10),  [vL  »0,  so  that  at  time  t0 
the  jump  conditions  (4.20c)  and  (4.20d)  reduce  to  the  corresponding 
conditions  for  the  equilibrium  problem  as  given  by  Erlcksen  [1], 

Although  there  are  local  existence  and  uniqueness  theorems  for 
nonsingular  differential -del ay  systems  and  for  singular  systems  with 
known  delays  [22],  these  do  not  apply  to  (4.20).  We  shall  therefore 
take  for  granted  the  existence  of  a  solution  to  our  system.  Moreover, 
we  shall  assume  that,  near  t*tQ,  the  unknowns  s  and.  f  are  asymptot 

P 

Ically  characterized  by 

s(t)n,J  Mt-t/*  ,  (5.11 

k*0  K  p 

f(,)' V|/k<-z-\fte)"l‘  .  (5-12 


where  the  constants  s, ,  m,  f,,  and  n,.,  which  are  to  be  determined,  must 
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be  such  that 


m  >  1 ,  n.  >  0,  < 

o  10 


°.  f0ahyji  ♦(B1)-B2),  (5.13) 


m j  > ,  nj>nk  for  j>k  . 


(5.14) 


We  shall  determine  m  .  n. ,  s^  and  f.  explicitly,  so  that  the  dominant 

0  10  I 

terms  for  t-*-ta  in  (5.11)  and  (5.12)  will  be  known. 

D 

For  simplicity,  we  restrict  attention  to  loadings  for  which 


and  materials  such  that 


oo(t8)>0  , 


W**  (8-j )  <  0  ; 


(5.15) 


(5.16) 


the  latter  of  these  conditions  ensures  that  the  stress-strain  curve  is 
locally  parabolic  (rather  than,  say,  quartic)  near  e  »  S-j  .  We  have 

a„(t)  »aB  +  a0(tg)(t- tfl) +0[(t-ta)2]  as  t-»-tft  . 


(5.17) 


Then  by  (3.9)  and  (4.19),  the  discontinuous  function  eQ(t)  is  given 


near  t  *  tg  by 


r 


eo(t)  *  ^ 


51  ■  £ W" ('ffjf  j  (tg-t)*  +o[(tg- t)*]  .as  t-tg  , 


(5.18) 


6,+ 


a0(tB) 


(t- tj  +  o  (t  -  t_) 


as  t-*-t. 


3* 


i  ■  -j 


t-"*r 

r..  ^A 


^See  (5.10) 
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For  use  In  (4.20e),  we  require  the  appropriate  approximation  to 
♦(e)  for  e-8"  .  From  (2.18)  and  (2.19), 


♦(e)  *  ♦( 8-j ) 


/p  ’ 


(5.19) 


so  that  making  use  of  (2.15)  and  (5.16)  we  obtain 


♦(e) 


♦(0,) 


2 

T 


^-W"(8,)/P 


( B-j  -e) 


o[(B, 


e)4] 


as 


(5.20) 


Since  *r(s(t),t)  occurs  as  an  argument  In  (4.20),  we  require  Its 

behavior  for  t-*t”  .  To  this  end,  one  first  observes  that,  by  (4.20h), 

0 

(2.15),  (5.16)  and  (5.18)^ ,  near  t  ■  t*  * 

s(t)'v-  ^(tg- x)^(t- t)  ,  (5.21) 


where  the  constant  ICj  Is  given  by 


-■[ 


-2lT(B1)o0(tB) 


]'• 


(5.22) 


A  detailed  analysis,  making  use  of  (5.11),  (5.21)  and  (5.22),  shows  that 
the  dominant  behavior  of  x(s(t) ,t)  depends  on  the  value  of  the  unknown 
exponent  rn  in  (5.11)  as  follows:  as  t-*-t*  , 

0  P 
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•  s  5/4  4in  /5 


4m./ 5 


Vhcf r  <*“V  °  +oC(t"V  °  3  If  V  5/4  , 


x(s(t),t)  •< 


tg-X  (t- tg)  +o(t- tg) 


if  mQ  *  5/4 


-s 


4m. -4 


V 


a  fH.-'T  4m  -4 

V(Tcf}  +°C(t  -  tg)  0  ]  if  mQ>  5/4  . 


Here  the  constant  X  is  the  unique  positive  root  of 

x5/4  +  x1/4+  *o.0. 

*1 


For  convenience,  we  define  constants  q*q(mQ)  and  K2  ■ 


4m  /5  for  m  <  5/4  , 

o  o 


q  ■  ^  1  for  mQ  ■  5/4  , 


^  4m0-4  for  mQ>  5/4  , 


and 


f(-s  /K])5/4  for  mQ  <  5/4  , 


*  < 


for  m  *  5/4  , 

o 


for  m  >  5/4 
o 


Then  (5.23)  can  be  written  as 


tg- T*r:K^(t  -  tg)q +o[(t  -  tg)q]  as  t-*-tg 


(5.23) 

(5.24) 

by 

(5.25) 

(5.26) 

(5.27) 


Since  by  (5.13)  we  have  mQ  > 1 ,  it  follows  that 


We  are  now  able  to  determine  the  behavior  near  t = tg  of  the  functions 
W'  ( eQ ( t ) ) ,  £0(t),  d> ( eQ ( t ) )  and  s(t))  which  occur  in  (4.20): 

W'U0(t))  - cr0 ( t )  *o6-o0(t^K2(t-tB)q  +  o[(t-t  )q],  (5.29] 


£0(t)*l-[-250(t8)K2/r'(B1)]J(t-t6)‘,/2  +  o[(t-t6)q/2]  ,  (5.30) 

♦(e0(T»  -  *(8, )  -  2/3(-W“  (61))r>p-i(2io(te)t2)i(t  -  tg)3q/4 

+  o[(t-tg)3q/4]  ,  (5.31) 

•  .  . 

*62+  "V’6,  (t~  V  +  -yjp  SoW^‘ V  ° 

+  o[(t-t6)m°]+OC(t-tB)2]  .  (5.32) 

Furthermore,  (5.11),  (5.12)  together  with  the  binomial  expansion  then 
yield 

f(s(t) -^/f  t>+  f(-s(t)  -^Tt>  =  2f0  +  2f,(>/|)"1(t- v"1 

♦  2f„(./l )"2  <t  -  t„)"2  ♦  of(t  -  tj"1"!"!4",-1  *n2^1 


(5.33 
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n,-l  nn+ino-1 


-f(s(t)  +f  (-s(t)  -^|t)  "2f1n1s0(^|)ni  (t-t6)  1 

rr  n.-l  n.+m,-l  rr  n?-l  m  +n?-l 

^WlWf*  (t“V  +  2f2n2soW {t“t6) 

,2mo+n^-2,mo+n2-l| 


[min{n,+ffl,-l  ,2m  +n,-2,m  +n,,-l 1 

(t-t6)  11  01  0  2  J 


)  ( 


+  0|  (t- 16) 


m1n{3tn0+n^-3,m^+n2-f  »m0+n3“U 


i 


(5.34) 

Equations  (5.33)  and  5.34)  make  it  possible  to  determine  the  behavior 
near  t = tQ  of  the  various  physical  quantities  which  suffer  jumps  across 

P 

the  phase  boundary.  Thus  from  (4.20e),(5.33) ,  (5.32)  and  (5.31) 


fr  n,+l  n.  rrn?+1  n 

W-Wf)  1  *2W|)  f2(t-V 


+  -2-5-  U-t6)4<-^(Bl>>Vi<2i0<tB)K2Wt-tB)3q/4 


/p  E 


+  o 


min(n.|-Hno-l  ,n2,l  ,3q/4 3 


(5.35) 


while  (4. 20f ) ,  (5.34),  (5.32).  (5.30)  and  (5.11)  give 


m  -1 


+0 


in  i  m* 

[e]s  =  (82-B1)moso(t-tg)  0  +(02  -  B1)m1SQ(t  -  tg) 

tvoK't8>vw"'(6i))J(t-tB),/2tm°'1 

min{n^+2mo-2,mo3 


m,  -  1 


m1n(n,,+m  -1 ,  q/2  +  m  -1  ,m, -1 } 
20 


(5.36) 


From  (4 .20 g )  ,(5.34) ,  (5.32)  and  (5.29)  we  have 


,  rr  n .+1  n;+m  -1  [7  n.+l 

-  [W7  ]  *  2f1n1sQ(^-)  ’  tg)  +2flnlsl  p)  (t-t 


n.+l  n,+  m,-l 

I  1  (t-tB)  1  1 


rr  n,+l  m  +n^-l  a  (tfl)  *„&«(*«) 

+  2f2"2so(>/|>  «*-V  -apJ-(t*te,*-^J-<t-ta> 

_  °0(t6>K2  (t  _  tg)q  +  0  r  _  -"o-U  j 


+  0  (t-tg) 


Finally  (5.35)  and  (5.11)  give 


m1n{3m0+n|-3,in^+n2-l  *m0+n3’U 


n.+m  -1 
,)  1  0  +  2m 


(7  n.+l  n,+rn  -1  17“  n,+l  n. 

wi-v.  s/|>  V‘-V  tVi(4)  fi(t-V 

[T  n.,+1  n9+m  -1  vA^a)  mA 

+2Vo(>/F>  ¥‘-V  +^r-(t-v 

9  xx  3  -x^-  +  m-1 

["  nilntHj-Hnj-l -1  ,n1+2m()-2,m0,T^'Hn0-l] 

+  0  I  (t  -  tg) 


n^-Hn^l 


+  o  (t-tg) 


It  Is  now  necessary  to  enter  the  two  shock  conditions  (4.20c) 
and  (4.20d)  and  Investigate  the  possible  dominant  balances  among  the 
various  terms  as  t  +  tg.  This  process  makes  use  of  the  asymptotic 
results  given  above  and  requires  a  detailed  analysis,  too  laborious  to  be 
included  here,  of  various  possible  cases.  One  Is  thus  led  to  the  values 


mQ  *  3/2  ,  n1  -  1/2 


(5.39) 


for  the  leading  exponents  and  to  the  formulas 


s 


o 


3<B2-6i )i0(t6 


1* 

•  > 


(5.40) 


Thus  the  dominant  terms  in  the  expansions  (5.11)  and  (5.12)  are  determined: 

* 

s(t)  ^sQ(t  -  tg)  ,  (5.41) 

ts)S  ,  (5.42) 

where  sQ  and  f^  are  given  by  (5.40).  Although  we  have  not  pursued 
In  full  detail  the  higher  order  corrections  to  (5.41)  and  (5.42),  we 
conjucture  that  the  values 


mk  ‘  7T  *  nk  "  7 


(5.43) 


for  the  exponents  In  (5.11)  and  (5.12)  provide  for  a  consistent  expan¬ 
sion  scheme. 


6.  Large-Time  Behavior  of  the  Phase  Boundary 

We  now  undertake  an  asymptotic  analysis  to  determine  the  large¬ 
time  behavior  of  the  propagating  phase  boundary.  The  loading  condition 
( t )  >. 0  suggests  that  we  seek  solutions  such  that 


s(t)  <0  for  t  >  tfl 


(6.1) 


The  validity  of  (6.1)  In  some  neighborhood  of  time  tfl  follows  from  the 

P 

results  of  the  previous  section.  Although  we  conjecture  that  (6.1)  fol¬ 
lows  from  the  governing  equations,  either  (4.20)  or  (4.32),  we  have  been 
unable  to  confirm  this.  We  shall  show,  however,  that  assumption  (6.1) 
leads  to  physically  reasonable  behavior  of  the  solution  to  our  problem. 
From  (4.15),  (6.1),  and  the  short-time  analysis  we  have 


-£■< 


s(t)  <  0  ,  t  >  tfl 


(6.2) 


Moreover,  since  s(t)  is  monotonic  decreasing  and  bounded  below, -it  fol¬ 
lows  that  It  approaches  some  limit  as  t  -*■  od  .  Thus  we  define  the 


)hase  boundary  speed  a  by 
a  *  -11m  s(t) 

t  -*■  OD 


hence 


as/F  • 


(6.3) 


(6.4) 


The  acoustic  speed  c(x,t)  Is  constant  on  each  negatively-sloped 
characteristic  ray  In  moreover,  c ( 0 , t^ )  =  0.  Thus  (6,1)  guarantees 
that,  with  Increasing  time,  the  phase  boundary  will  intersect  characteristic 
rays  In  which  arise  on  x=0  at  progressively  earlier  times;  see  Fig.  8. 
This  geometric  state  of  affairs  is  described  by  the  inequalities: 


.<  »>  v 


— »  'Jm  l.  m  .  »  -  »  .  >  „ 


.  *.*  *,„*»•  IJVa. 


O  V*  •.*  %' 
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^  r(s(t),t)  <  0  ,  t  >  t$  ,  (6.5) 

and 

c(e0(x(s(t).t)))  <  -s(t)  ,  t  >  tg  .  (6.6) 

The  validity  and  equivalence  of  (6.5)  and  (6.6)  also  follow  from 
(6.1),  (4.20h)  and  the  short-time  analysis  of  the  previous  section.  From 
(2.15),  (2.18),  (2.19),  (2.32),  (2.36),  and  (6.5)  we  infer 

^  eQ(r(s(t),t))  <.  0  and  ^  <J>(eQ(T(s(t),t)))  <  0  ,  (6.7) 

so  that  e(s(t)",t)  and  v(s(t)",t)  are  decreasing  with  time.  We  have 
already  observed  that  the  first  inequality  in  (6.2)  guarantees  that,  with 
increasing  time,  the  phase  boundary  will  intersect  members  of  the  nega¬ 
tively-sloped  family  of  characteristic  rays  in  fc2  which 'Issue  progres¬ 
sively  later  on  x-0.1  This  alone  does  not  allow  us  to  draw  any  conclusion 
regarding  the  monotonicity  of  e(s(t)  ,t)  and  v(s(t)  ,t).  We  will, 
however,  be  able  to  establish  such  a  monotonicity  result  once  we  deter¬ 
mine  some  of  the  properties  of  the  unknown  function  f(z).  Meanwhile, 
we  consolidate  (6.2)  and  (6.6)  into 

->/£<  s(t)  <  -c(eQ(T(s(t),t)))  ,  t  >  tg  .  (6.8) 

Thus  the  speed  of  the  phase  boundary  Is  confined  to  the  interval  between  the 
different  acoustic  speeds  on  each  side  of  this  boundary.  We  note  that 
(6.8)  is  precisely  the  entropy  condition  given  by  Lax  [18], 

It  is  apparent  that  as  t  -*■  od  one  of  the  following  two  alternatives 
must  occur. 

^See  the  discussion  following  (4.18). 
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Altemative  I:  The  phase  boundary  In  the  x-t  plane  does  not  cross 
the  Initial  characteristic  x*-c(0)t  and  so  is  confined  to  Rg. 

For  this  to  occur  it  is  necessary  that 

-c(0)t  <  s(t)  for  t  >  tg  .  (6.9) 

Since  by  (6.3),  s(t)  is  asymptotically  linear,  (6.9)  will  hold  if 
and  only  if 

a  <  c(0)  =  .  (6.10) 

Thus  for  large-time,  the  curve  x*s(t)  will  be  asymptotically 
parallel  to  a  limiting  characteristic  C_  ray  in  Rj.  We  shall  let 
denote  the  time  at  which  this  limiting  characteristic  arises 


on  x*0.  Hence, 

t  -  11m  r(s(t),t)  , 

t  -*•  QD 

(6-11) 

and 

“  "  C<e0<T,B» 

(6.12) 

Alternative  II:  The  phase  boundary  in  the  x-t  plane  crosses  the  ini¬ 
tial  characteristic  at  some  time  t^  which  satisfies  the  equation 

s(tL) =  -c(0)tL  .  (6.13) 

Then 

s(t)  <  -c(0)t  for  t  >  tL  ,  (6.14) 

so  that  after  time  tL  the  phase  boundary  is  the  leading  disturb¬ 
ance.  From  (6.14)  it  follows  that  this  alternative  will  occur  if 
and  only  if 
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ct  >  c(0)  =  .  (6.15) 

It  Is  clear  from  (6.10)  and  (6.15)  that  whether  or  not  the  phase 
boundary  becomes  the  leading  disturbance  depends  on  the  asymptotic  phase 
boundary  speed  a.  Thus  a  determines  which  of  the  foregoing  alternatives  occurs 
To  ascertain  the  large- time  behavior  of  the  phase  boundary,  we 
shall  study  the  large-time  behavior  of  solutions  to  the  pair  of  differ¬ 
ential-delay  equations  (4.22).  A  major  difficulty  In  any  such  asymptotic 
scheme  stemsfrom  the  lack  of  an  explicit  formula  for  -r(s(t),t)  when  the  phase 
boundary  Is  In  V  We  faced  the  same  difficulty  In  connection  with  the 
short-time  behavior  where  an  asymptotic  analysis  led  to  (5.23).  For 
large  time,  another  approximation  scheme  Inmedlately  suggests  Itself.  In 

this  case  we  wish  to  consider  t'» t  >x(s(t),t),  so  that  from  (4.20h), 

B 

c(c0(t))  .  (6.16) 

Now  with  (3.2),  the  discussion  In  section  2  shows  that  c(eq(t))  Is  con¬ 
tinuous  and  monotonlcally  decreasing  as  a  function  of  x.  If  c(e0(x)) 

Is  strictly  decreasing,  then  It  Is  Invertible.  By  (2.34),  this Invertlbll- 
Ity  Is  assured  If  both 


aQ(t)  >  0 

0  <  t  i  ts 

(6.17) 

w"'(e)  <  0  .  , 

o  <  e  <  61  . 

(6.18) 

Note  that  (6.18)  tightens  assumption  (2.38)  by  requiring  that  the  mate¬ 
rial  Is  not  only  "non-hardening  In  the  first  phase"  but  also  genuinely 
nonlinear  or  "softening  In  the  first  phase."  In  the  remaining  analysis 
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we  shall  assume  that  (6.17)  and  (6.18)  hold.  Thus  we  may  define  an 
Inverse  ft  for  c(eQ(*))»  where 


and 


n:  [0»c(0)]  [o.tB]  , 
z  *  c(e  (fl(z))) 


t  *  ft(c(e0(t)) 


By  (6.16)  It  now  follows  that 


Oiz<  c(0)  t 
0  <  t  <  tg 


r  (6.i9) 


J 


for  t  »  t. 


provided 


0  <  llM  <  c(0) 


(6.20) 


(6.21) 


If  the  phase  boundary  becomes  the  leading  disturbance,  then 
Alternative  II  holds  and  s(t)  €  ftj.  for  t  >  tL;  consequently, 
t(s (t) »t)  *  0  for  t>tL.  Thus  a  large-time  asymptotic  study  of  (4.22) 
requires  an  explicit  approximation  for  x(s(t),t)  only  for  Alternative  I, 
the  alternative  where  the  phase  boundary  does  not  become  the  leading 
disturbance.  Then  (6.1),  (6.3),  and  (6.10)  yield 


0<-s(t)<c(0)  ,  t  >  tg  .  (6.22) 

Since  s(t)  Is  asymptotically  linear,  (6.22)  Implies  (6.21)  for 
large  times.  Thus  when  Alternative  I  holds  we  may  appeal  to  (6.20)  In 
the  asymptotic  analysis.  We  note  that  for  such  an  analysis,  0 
need  only  be  defined  on  some  one-sided  neighborhood  of  a,  say  (a -6, a]. 
Hence  we  could  weaken  (6.17)  and  (6.18)  to 
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W  >  0 


•*' <<=„<Too»<0  • 


(6.23) 


However,  as  t  Is  not  known  In  advance.  It  is  convenient  to  retain  (6.17) 
oo 

and  (6.18). 

When  Alternative  I  holds,  (6.20)  furnishes  the  leading  term  In  an 
asymptotic  expansion  of  x(s(t),t)  for  large  t.  This  can  be  seen  by 
considering  the  following  Iteration  scheme  suggested  by  (4.20h): 


rn+l  *  • 


T0*° 


(6.24) 


for  which  x(s(t),t)  Is  a  fixed  point  and  (6.20)  Is  the  first  iterate. 
Equation  (6.24)  generates  a  formal  series  solution  to  (4.20h)  of  the 


r(s(t),t)  »  l  Ak(^^-)  t'k 


(6.25) 


k  / 

where  A.(z)  Is  a  finite  sum  of  products  of  z,  n(z),  4r  n(z)*  . 

dzk 

The  first  three  A^'s  are  given  by 
Aq(z)  =  n(z)  , 

A,(z)  *  2 n(2)  . 


V 2)  ■  22I!(2)(^a^52  -  2(n(2))Z 

+  J  Z2(n(z))2  . 

£  dz^ 

Differentiation  of  the  second  of  (6.19)  yields 


(6.26) 


rrW 


•V-V 

m 


■  -.i 

•v 

:v-c; 


£n(z) 


c(e0(n(z)))  e0(n(z)) 


(6.27) 


1  V  *  *•.  -\  •*. 

"oV'*.  •’/v  V  '.' •• ' 
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It  then  follows  from  (2.18)  and  (2.32)  that  (6.17)  and  (6.18)  are  pre¬ 
cisely  the  conditions  necessary  for  dfl/dz  to  be  finite.  It  similarly 
follows  that  all  the  A^z)  are  finite,  although  not  necessarily  uniform¬ 
ly  bounded,  so  that  (6.25)  Is  Indeed  an  asymptotic  series. 


From  (6.17)  and  (6.18)  It  follows  that  eQ(t)  and  c it) 

are  Invertible  In  the  first  phase;  their  Inverses,  denoted  by 

e-1  and  •  c"\  are  then  defined  on  0  <  e  <  g,  and  0  <.  c  <.  c(0),  respec- 
o  i 

tlvely.  We  then  have 


fl(z)  *  e^1(c"1(z)) 


,  (6.28) 
and  hence  arrive  at  the  following  asymptotic  formulas  for  Alternative!: 


T(s(t),t)  -  +  0(1) 

e(s(t)\t)-e#(T(s(t).t))  ^  C_1(^1H)  . 

M'(e(s(t)‘,t))  -  H'(e0(T(s(t).t)))  -v  u'(c''(^M))  , 

v(s(t)’,t)  «  *(e  (r(s(t).t)))  *  . 


(6.29) 

(6.30) 

(6.31) 

(6.32) 


Entering  (4.22)  and  (4.23)  with  (6 .29)-(6 .32) ,  we  arrive  at  the  asymptotic 
governing  equations  when  the  phase  boundary  does  not  become  the  leading 
disturbance  (Alternative  I): 


tqeWa 


S-*:l 


£3* 


I-  •'"’I 
---VO 

m 

V-V- 


>  dM\lb  ojt+j*  s(t)>  - 


'p'  o* 


(6.33a) 


(-i(t)^f)  t)  c-1!^) 


(6.33b) 


Conversely  entering  (4.22)  and  (4.23)  with  eQ(T(s(t) ,t))  *  0t  we  arrive 
at  the  exact  governing  equations  when  the  phase  boundary  becomes  the 
leading  disturbance  and  t  >  tL  (Alternative  II): 

>/F  (*(t)  %/f)  «-*<»>  \Jf< n ♦ 7=  <*(»)  +i/f>  «.(t+ji*(t)) 


D 


<i<t>+WF)  ■  °  • 


and 


JF  s/F>  ->/!*)  -  *  f  -  °  • 


(6.34a) 


(6.34b) 


In  proceeding  with  the  analysis  of  (6.33),  (6.34)  we  shall  treat 
four  separate  cases: 


(I) 

W"(0UE  , 

a  (t)  -*■  +cd 

0 

9 

(6.35) 

(II) 

w"(0)>E  , 

a_(t)  -*■  a_  < 

O  00 

+00  , 

(6.36) 

(III) 

w"(0)<E  , 

a_  ( t )  -*■  a  < 
0  OD 

+a>  , 

(6.37) 

(IV) 

W"(0)<E  , 

aQ(t)  -*■  +ao 

• 

(6.38) 

Since  the  end  condition  (6.17)  implies  that  lim  a  At)  <. +oo  ,  the  four  cases 

t+oo  0 


T 
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llsted  above  are  exhaustive.  Although  It  seems  likely  that  W//(0)  >.  E 
for  real  materials,  we  shall  treat  cases  (III)  and  (IV)  as  well  for  the 
sake  of  completeness.  We  note  from  (6.4)  and  (6.15)  that  the  condition 

W"(0)  <  E  (6.39) 

Is  necessary  for  the  phase  boundary  to  become  the  leading  disturbance. 

We  have  so  far  remained  silent  on  the  behavior  of  f(z),  the  re¬ 
maining  unknown  In  the  differential -delay  equations.  We  first  show  that 
f(z)  Is  monotonlcally  decreasing  on  Its  Interval  of  definition 

(-0>,-J£  tfl].  For  this  purpose,  as  well  as  for  the  forthcoming 
v  p  p 

analysis,  we  will  need  the  following  lemma. 

Lemma:  If  W(e)  satisfies  (2.14),  (2.15),  (2.16),  (2.17),  (4.2),  and 


(6.18),  then 

M'(c*’(E))  -  Ec'\o  -  D>  0 

(6.40) 

i 

|  where 

1  0  <  5  <  m1n{^— ^  ,  yfi)  . 

(6.41) 

We  note  that  m1n{  yj  }  Is  yj  In  cases  III, 

1  In  Cases  I,  II. 

IV;  It  Is 

Proof  of  Lemma: 

Let 

H(S)  =  W/(c”1(E))  -  Ec_1(e)-D  , 

(6.42) 

then  since  c'1(0)  =  81  It  follows  that 

H(0)  =  W/(8] )  -  ES-j  -  D  >  W/(61)-  E62-D  *  0  . 

(6.43) 

Since  H(0)  >  0,  we  need  only  demonstrate  that  H(0  Is  monotonlcally 


-9} 


increasing.  Writing 


h'(e)  -  (w'V'un-  e)  , 

we  examine  each  factor  of  H* (£)  Individually.  Now 

J _ 


d  c_1(£) 


3T 


de 


c(e) 


e  *  c"1 (^) 


As  for  the  first  factor  of  (6.44), 


WM(0)  -E  <0 
E-  E  -  0 


and 

^  (Wm(c"1(E))  -  E)  -  w"'(c-'(5)»  ^  c*’(E)  > 
Now  (6.46)  and  (6.47)  give 

W,,(c’1(0)-  E  <  0  . 


Hence  we  obtain  the  monotonicity  result 


H*(E)  >  0 


which  completes  the  proof  of  the  Lemma. 


Corollary:  For  W(e)  as  in  the  above  Lemma  and 


s  0  . 


for  W*  (0)  <  E 
for  W*(0)>E 


c*  defined  by 
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*  J  0  for  WM(0)  s  E 

l  C_1(>/J)  for  w"(°)i  E 

W'(e)  -  Ee  -  D  >  0 


(6.50) 


(6.51) 


when 


e*  <  e  <  61 


(6.52) 


We  now  show  that  f'(z)<0.  Differentiating  (4.22b)  and  using 
(4.23b)  we  find  after  some  algebra  that 


1 


2VfwH{t))i 


litl 


p( 


JF-* 


S(t)) 


[- 


(eQ(x(s(t),t)j 


-  EeQ  (x(s(t).t))  -  D 


c(e  (x(s(t),t))))  (s(t)  +  c(e  (x(s(t) 


*  ^  e0(T(s(t) 


,t))  •  . 


] 


(6.53) 


Later  on  we  shall  also  require  f '(-s(t)  - ^=-t),  which  we  display  at  this 
juncture.  From  (4.22a),  (4.23a),  we  obtain  after  elementary  manipulations. 


f (-s(t) 


f(z) 
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z  =  -s(t) 


1 


lit} 


2 


-#• 

- [w'(e  (r(s(t).t))) 

*  »  -  \  \  ■ 


P  P 


-  Ee0(x(s(t)tt)) 


►  ,t))  -  dJ 

[(-^f +c(£0(T(s(t),t))))  (s(t)  +c(e0(T(s(t),t))))J^-e0(x(s(t),t)) 


+  2(s(t)+>/i)2  i0(t  +  ^|s(t))j 

Now  recall  that 

ieo(T(s(t),t))  i° 


(6.54) 


(6.55) 


In  addition,  if  W*(0)  >.  E,  the  phase  boundary  cannot  become  the  leading 
disturbance;  accordingly,  (6.7),  (6.30),  (6.3),  (6.45),  and  (6.4)  yield 


6<|  >  e0(T(s(t),t))  >  c_1(a)  >  c_1(^)  when  w"(0)>  E.  (6.56) 

Consequently  e0(T(s(t),t))  always  occurs  in  the  interval  given  by  (6.52), 
so  that  the  Corollary  yields 

W#(e0(T(s(t),t)))  -Ee0(T(s(t),t))  -  D  >  0  .  (6.57) 

Examining  the  second  bracketed  quantity  in  (6.53)  we  have  with  the  aid 
of  (6.8)  that 


<4 


(J§-  c(e  (t(s ( t ) , t ) ) ))  (s(t)  +  c(e0(T(s(t),t)))) 


<  0,  t  >  tg .(6.58) 


Thus  applying  (6.1),  (6.8),  (6.57),  (6.58),  and  (6.7)  to  (6.53),  we 


conclude  that 


f'(s(t)  -y^t)  <0  ,  t  >  tg  . 

Moreover,  s(t)  <.  0  and  s(tg)  =  0  guarantees  that  t«-  s(t) 
Continuous,  one-to-one  mapping  of  the  interval  (-oo*  tg]  onto 
This  allows  us  to  conclude  from  (6.59)  that 


(6.59) 

•  Vf  ys  a 

(-00  »-yjp  tp]« 


f'(z)  <  0 


for 


-oo 


(6.60) 


thus  establishing  the  monotonicity  of  f.  It  Is  interesting  to  note  that 
it  Is  not  obvious  from  (6.54)  that  f;(-s(t)-  <  0,  although  of 

course  (6.60)  provides  this  result  as  well. 

Having  shown  that  f  Is  monotonically  decreasing,  we  now  demonstrate 
that  f(z)  is  bounded.  From  (5.10)  we  already  have  the  minimum  value  for 
f.  To  show  that  f  is  bounded  above,  we  appeal  to  (4.22b),  which  yields 


f(z)  < 


(t),s(t),t) 


< 


(6.61) 


In  view  of  (6.60),  (6.61),  lim  f(z)  exists,  and  so  we  may  define 

Z-*-  -CD 

f  s  11m  f(z)  .  (6.62) 

Z  -►-CO 

Combining  (5.10)  and  (6.62)  , 

l(^f*(8,)-62)  =  <  f(z)  <  ^  .  (6.63) 

The  preceding  discussion  concerning  the  behavior  of  s(t)  for  large 
positive  t  and  of  f(z)  for  large  negative  z  leads  us  to  seek  asymptotic 
solutions  of  (6.33)  or  (6.34)  in  the  form 


s(t)  =  -  at  +  q(t) 


(6.64) 


f<z>  *  f-a>  +  fl<z)  • 


(6.65) 


where 


q(t)  =  o(t)  ,  q(t)  *  o(l)  as  t  -*■  oo  , 


(6.66) 


f1  (z)  =  o(l )  as  z  -*■  -od 


(6.67) 


We  shall  give  formulas  determining  a  and  f  ,  and  shall  seek  more 


specific  Information  about  q(t)  and  f^(z). 

In  the  ensuing  analysis  of  (6.33)  or  (6.34)  we  are  faced  with 


determining  the  asymptotic  behavior  of  f(-s(t)  -  */^t)  ,  f(s(t) 


•4 


t) 


and  aQ(t+Jf  s(t))  as  t  oo .  Here  we  shall  use  the  fact  that  since 


f  and  aQ  are  monotonic,  their  asymptotic  behavior  as  their  argument  ap¬ 


proaches  some  value  (finite  or  infinite)  Is  simply  given  by  the  function 

1 


evaluated  at  Its  asymptotic  argument. 


Case  I:  W  (0)>.E,  qQ(t)  •»  oo 


Here  (6.33)  necessarily  applies  and  substitution  of  (6.64)  Into 
(6.33a)  yields  the  asymptotic  equation: 


4 

(-C  +  \ $  +  x\a)  -j(-c  ♦(«-’  M)  -  jpW/(c*’1  (c))  =  0 


(6.6ft) 

as  Is  the 


The  last  four  terms  are  bounded  for  all  a  with  0< 
first  term  ,  because  of  (6.63).  However,  if  a  second  term  is 


1 


This  procedure  was  tacitly  employed  In  obtaining  (6.30)-(6.32). 


unbounded.  Thus  we  conclude  that 


When  (6.69)  holds,  (4.17)  requires  that 


(6.1 


q(t)  >  yjf  tg  .  (6.: 

On  account  of  (6.69),  equation  (6.33a)  to  dominant  order  yields 
J|q(t)  f(-q(t))+-^q(t)  q(t))+^ 

•/pc 

c*1(>/F)-iw,(c'1(>if))  *°  *  (6-: 

while  (6.33b)  to  dominant  order  gives 

-V^(c',(A/I))+iw'(c",(>/F))  s0  •  (6-: 

From  (6.72)  and  (6.65)  we  obtain 


Turning  our  attention  to  (6.71)  we  note  that  the  last  three  terms  are 
0(1)  and  by  the  Lemma  do  not  sum  to  zero.  Since  q(t)  ao(l)  and  f  is 


bounded,  It  follows  that  to  effect  the  dominant  balance,  cQ(  ^/fqlt)) 


must  be  unbounded  as  t  -*■  co  .  This  requires  that 


11m  q(t)  =  +oo  , 
t  +  oo 


(6. 


ice  (6.70)  Is  satisfied  for  1 


and  that  the  leading  behavior  of  q(t)  Is  supplied  by  the  solution  of 
q(t)o0 (^/fqlt))  >o  .  (6.; 

the  Inequality  being  a  consequence  of  the  Lemma.  Integration  of  the 
^bove  expression  yields 

j  a0<  ,/f  *)de  -  -  Ec',(  S  '  0]t  •  (6': 

Let  us  consider  as  examples  power- law  and  exponential  loadings: 

Example  1:  Asymptotic  power-law  loading 
Let 

aQ(t)  *  k-jtn  +  o(tn)  ,  k1  >  0,  0  <  n  <  oo .  (6.1 

Upon  substitution  Into  (6.76)  this  yields 

WfcwVl(>/f))-  Ec'1(>/I)' o]t 

'I  *  kl(t,n/2  ffFT  •  (6-: 

so  that 

t(t)  ^  Ec'1(^f)-“))n+'  ^  ■  <6‘: 


Example  2:  Asymptotic  exponential  loading 
Let 


k-t  k,t 

o^(t)  =  kg  e  + o(e  )  »  kg  >  0,  k^  >  0  . 


(6.1 


Substitution  Into  (6.76)  now  yields 


^WV’^I-EC-1^-03* 


(6.81) 


whereupon 

q(t)  *  Y^yff  in  1  ‘ 


(6.82) 


From  these  examples,  we  see  that  the  faster  oQ(t)  tends  to  Infinity 
the  smaller  Is  the  dominant  behavior  of  q(t).  Hence  the  phase  boundary 
approaches  Its  asymptotic  speed ^  faster  for  faster  loading  rates. 

Since  the  final  round  of  the  above  analysis  Involved  a  dominant 
balance  among  the  o(l)  terms  In  (4.22),  further  corrections  to  f(z)  and 
q(t)  will  involve  consideration  of  higher  order  terms  from  (6.30)-(6.32) 
Because  this  requires  considering  the  oO^)  term  In  (6.25),  It  appears 
that  higher  order  analysis  would  become  an  order  of  magnitude  more 
tedious. 

Case  II:  ldw(0)>  E  ,  o0(t)  -*'O0D<a > 


Again  the  phase  boundary  cannot  become  the  leading  disturbance  so 
that  (6.33)  apply.  We  can  dismiss  the  possibility  a  =  »  since  In 

that  case  the  first  two  summands  In  (6.33a)  are  o(l)  and  we  are  left 
with  the  0(1)  terms  which  must  satisfy 

^  (D+Ec'’(^|)  -  Me ''(-v/f)))  -  0  •  («•«) 


contradicting  the  Lemma.  Hence  we  conclude  that 
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0  <  a  <  ^  .  (6.84) 

Consequently,  not  only  Is 

f(s(t)  _‘^t)  *  f-oo  as  (6.85) 

but  also 

f(-s(t)-^|t)  ^f((a-^|)t)  ^  f-aj  as  t-«.  (6.86) 

Entering  with  (6.85),  (6.86),  and  (6.64)  Into  (6.33),  we  are  lead  to  the 
asymptotic  equations: 


and 


(-*  +  Ji)  f  m+  —  ( -a  +  /^)ct  -  —(-a  + 

VP  VP  -00  \P  00  2\P 

+  ? “a/I  c‘,(a)  '  I(-°  (*) )  ‘  £  WV1  (a))'0,  (6.87) 

>/f  f-a>'  ^  +  aC"’(a) 

-  jfo  +  Jf)  +(c‘'(a))  +  ^  H’(c_1(ot))  •  0  .  (6.88) 


Eliminating  f  between  (6.87)  and  (6.88)  we  arrive  at  an  equation  for  a: 
G(a)  =  (a2--)  ^+lw,(c"1(a))-a2c"1(o)-4^  =  0  .  (6.89) 

p  t  p  t 

Any  root  of  G(a)»0  1nthe  Interval  0  <  a  furnishes  a  candidate  for 
the  asymptotic  phase  boundary  speed.  We  now  show  that  there  is  a  unique 
root  a  of  (6.89)  in  the  above  Interval.  To  this  end  note  that 


s(°>  *  p  a®  P  w  (c  '<°»  s-F(cToo"ao^^  <  0  • 


(6.90) 


6(^f)  "  P  (W'(c’1(>/|))-  E  c"1(>/f)-  D)  >  0  ,  (6.91) 

where  the  Inequality  In  (6.91)  follows  from  the  Lemma.  Thus  we  are 
guaranteed  at  least  one  root  In  the  desired  Interval.  To  establish  the 
uniqueness  of  this  root,  we  demonstrate  that  6(a)  Is  strictly  Increasing  on  the 
Interval  0  <  a  <  Ti  .  From  (6.89), 

§  '  T  (aa>-  Ec',(a>  *  D>  +  <?  «  V’(«»  -  «2)  |f  «•’(<.)  .  (6.92) 


Whereas  (2.18)  gives 


(1hV'(o))-  a2)  ■  0  , 


which  implies  that  (6.92)  may  be  rewritten  as 


(6.93) 


a!  *  T  (fra’V*  CaB  -  W'(c-'(a))]+  [wV(a))  -  Ec"'(a)  -  D])  .  (6.94) 


'/--l 


'/  -1 


Since  each  of  the  bracketed  quantities  In  (6.94)  Is  positive  on 
0<a<  ,  It  follows  that 

f  >  0  for  0<o<^. 


(6.95) 


Consequently  (6.89)  furnishes  a  unique  value  for  the  asymptotic  phase 
boundary  speed  in  the  interval  0  <  a  <  0ne  may  now  find 

from  either  (6.87)  or  (6.88). 


From  (6.89)  It  follows  that  the  larger  ,  the  larger  Is  the 
asymptotic  phase  boundary  speed.  Moreover,  as  the  last  three  terms 
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of  (6.89)  are  bounded,  the  unique  root  approaches  as  tends  to 
Infinity.  This,  of  course,  Is  to  be  expected  from  the  results  for  Case  I. 

In  contrast  to  Case  I,  we  have  obtained  here,  from  the  0(1)  terms 
In  (6.33),  only  the  leading  asymptotic  terms  for  s(t)  and  f(z)  which  are 
a  and  f  .  Hence  the  first  corrections  to  both  s(t)  and  f(z)  will 
undoubtedly  require  consideration  of  the  o(-^)  term  In  (6.25). 

We  now  consider  materials  Which  admit  the  possibility  that  the 
propagating  phase  boundary  becomes  the  leading  disturbance,  that  Is, 
materials  for  which 

W"(0)<  E  ,  c(0)  * y[ <  yjf  •  (6.96) 

It  then  follows  from  (2.14),  (2.15),  (2.17),  (4.2),  (6.18),  and  an  ap¬ 
plication  of  the  mean  value  theorem  that 

D  <  0  .  (6.97) 

This  Inequality  Is  also  apparent  from  the  graph  of  an  appropriate 
stress-strain  curve  that  obeys  (6.96). 

Since  the  phase  boundary  will  not  become  the  leading  disturbance 
If  a  <  c(0)  =  (see  (6.10)),  (6.33)  applies  when 

0  <  a  £  c(0)  .  (6.98) 


D=  W'teg)  -  EB2  <  W'(B1 )  -  EB,  =  W*(0)  +  W'^O^  +  \  W,,#(s)B*  -  EB-j 
=  (Ww(0)  -E)B1  +  ^  W#,,(F)Bi<(w',(0)  -  E)  B-,  <  0,  where  0  <  e  <  Bj  . 


1 
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Conversely,  the  phase  boundary  becomes  the  leading  disturbance  if 
a  >  c(0), so  that  (6.34)  applies  when 


c( 0)  <  a  £ 


(6.99) 


Case  III:  W*(0)  <  E  ,  oQ(t)  -*■  oQD<  oo 


Let  us  suppose,  prior  to  the  analysis,  that  the  phase  boundary 
does  not  become  the  leading  disturbance.  Then  we  can  once  again  appeal 
to  the  analysis  which  yielded  (6.89),  provided  we  also  satisfy  (6.98). 

The  roots  of  G(a)  *  0,w1th  a  satisfying  (6.98)  furnish  candidates  for 
the  asymptotic  phase  boundary  speed.  Once  again, 

<3(0)  <  0  ,  ^|>0  for  0  <  o  <  c(0).  (6.100) 

Hence  (6.89)  has  a  unique  root  In  the  interval  given  by  (6.98)  if  and 
only  If 

o  2 

G(c(0))  *  (c(0)2-|)  -f- -  >  0  .  (6.101) 

Therefore,  a  necessary  condition  for  the  phase  boundary  not  to  become  the 
leading  disturbance  is  that 


GO 


<  -D-ejOIl 

f-c«» 


(6.102) 


Now  let  us  suppose  first  that  the  phase  boundary  becomes 
the  leading  disturbance.  Again  we  denote  the  time  at  which  this 
occurs  by  t^.  Then  (6.34)  applies  when  t  >  t^.  Here  the  asymptotic 


vw  7  V  •>  .  -  vvv»  /•  *. 

t  *  •*  •  O  •*  •  .4  ^  ,  *  i  *  *  ,  4  ,  »  ,  <  •  *  .*■,  »  .  i 
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7p  =  0  * 


(6.103) 


contradicting  (6.97).  Hence,  c(0)  <  a  <  ,  so  that  we  again  have 

(6.85)  and  (6.86).  The  asymptotic  equations  provided  by  (6.34)  are 

>/!<-<■♦#  +  ■°-  (6-104) 


>/§<■♦>/§> 


f  ■ !  o  =  °  . 

p  “00  c  p 


(6.105) 


Elimination  of  f  m  between  (6.104)  and  (6.105)  leads  to  the  equation 
-0D 


(-a2+|)  a2D  *  0  , 


(6.106) 


which  has  the  single  positive  root 


CD 


r 7i  * 


(6.107) 


(6.97),  a  <  yj^  .  Finally,  the  first  Inequality  In  (6.99)  requires 


(1  -JT-)' 

CD 


>  C(0)  , 


(6.108) 


phase  boundary  speed  must  lie  In  the  Interval  given  by  (6.99).  As  before 
we  dismiss  the  possibility  that  a  ■  since  then  the  asymptotic 
equation  resulting  from  (6.34)  Requires  that 


or  equivalently. 
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a_  > 
CD 


-D  c(0) 

f-c(0): 


(6.109) 


Inequalities  (6.102)  and  (6.109)  reveal  the  circumstances  under 
which  the  phase  boundary  becomes  the  leading  disturbance.  If  the 
ultimate  end  load  Is  less  than  or  equal  to  a  cut-off  load 


a, 


(6.110) 


the  phase  boundary  Is  confined  to  the  disturbed  region  and  the  asymp¬ 
totic  phase  boundary  speed  Is  given  by  the  unique  root  of  G(a)  =  0  In 
the  Interval  0  <  a  <.  c(0).  Conversely,  If  >oL,  the  phase  bound¬ 
ary  becomes  the  leading  disturbance  and  the  asymptotic  phase  boundary 
speed  is  given  by  (6.107).  We  note  that  an  application  of  the  mean- 
value  theorem  (In  a  manner  similar  to  that  which  led  to  (6.97))  yields 


al  »  °6  *  ao(V  • 


(6.111) 


If  we  define  a  as  a  function  of  via  (6.89)  when  (6.102)  holds, 
and  via  (6.107)  when  (6.109)  holds.  It  will  follow  that  a  is  a  continu¬ 
ous  and  increasing  function  on  <  aQ0<  00.  As  approaches  a^, 
one  has  a  -*■  0.  When  =  aL,  It  follows  that  a=  c(O),  so  that  in  the 
x-t  plane  the  phase  boundary  is  confined  to  the  disturbed  region,  but 
Is  asymptotically  parallel  to  the  leading  characteristic  curve 

x=  -c(0)t.  Finally  as  a  00  ,  a  approaches 

00  p 

Having  obtained  a,  we  find  f  when  <.  from  either  (6.87) 
or  (6.88).  Here,  as  in  Case  II,  further  corrections  to  s(t)  and  f(z) 


will  again  require  consideration  of  the  O(^)  term  in  (6.25). 

When  >  oL  we  obtain  f>a)from  either  (6.104)  or  (6.105).  The 
result  is 


-oo 


D 

2/pE 


1 


(<*+  \J-) 


s 


D 

2T 


0-^)1/2 

CD  .  « 


(6.112) 


p'  co 

Asymptotic  corrections  to  all  orders  for  f(z)  and  s(t)  are  found  *rom 
(6.34).  Since  the  analysis  is  not  difficult  and  the  result  is  interest¬ 
ing,  we  proceed. to  obtain  some  of  these  corrections.  Let 


a  (t)  *  0^+  a-j  (t)  ,  (t)  >  0  ,  o^  (t)  =  o(l )  as  t  -*■  oo . 

(6.113) 

Substituting  from  (6.113),  (6.64),  (6.65)  into  (6.34)  and  cancelling  the  pre 
vlously  balanced  o(l)-terms  eventually  yields  the  asymptotic  equations: 

(f— <D+  ~T  '  +  (“a  +  ^  V(a“  >/f) 4  ) 

+  ("a+>/|)  rai((1"  Vta)t)  =°  *  {6'114) 

-q(f)  ( <*+> /|)  f-| ( (-a -  yj^)t )  =  0  .  (6.115) 

Eliminating  q(t)  between  (6.114)  and  (6.115),  one  has 

(f-oo+  “T  ‘  f)(a+  fl((_a‘  VfW 

+  f-co(_a+>/|)  fi((a- y/fw 

al^“Vt° l)t)  =  0 


j 


+ 


» 


(6.116) 
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which  may  be  used  to  determine  the  asymptotic  behavior  of  f-j(z).  From 
(6.115)  the  asymptotic  behavior  of  q(t)  is  then  given  by 


.  •  (a  +  J|)  rr 

q(t)  ^  f - 8-  f 7  ( ( -ot  -  yjkt)  . 


(6.117) 


We  are  also  now  in  a  position  to  determine  the  asymptotic  behavior  of 
t 

I  f-|  ((-a  - -yj=-)e)ds,  where  t*  is  some  fixed  time  greater  than  t^.  If 
t* 

this  Integral  is  unbounded  as  t  -*■  oo ,  it  follows  upon  integration  of 
(6.117)  that  the  contribution  due  to  the  lower  limit  of  integration  t* 


Is  subdominant.  In  this  case  we  may  further  conclude  that 


(a  +  ^f)  r  If 

* — f — J  f -j ( ( -ct -  yjpa)d8  . 


-oo  ^ 


Consider  as  an  example. 


(6.118) 


oQ(t)  *  0^+  c^t""  +  o(fn)  ,  a1  <  0  ,  n  >  0  . 


(6.119) 


Then  (6.116)  suggests  the  asymptotic  form 


as  z  -*■  -oo  , 


(6.120) 


where  7-j  is  a  constant  to  be  determined.  Substitution  from  (6.119)  and 
(6.120)  into  (6.116)  yields 

”f-oo  /,  fo  \’n+1 

, _ ~^f(1'^a)  °i _  „ 

1  /7  -n+l  ItT  -n+l  cr  -  D  Hr  -n+1  <  * 

f-®((aWf>’  +(-a+V f)  )+JV(a+\f> 

(6.121) 

where  the  inequality  follows  from  a-j  <  0,  a  <  and  (6.112).  Further 


from  (6.117)  , 


7  rr  -n+1 

~  (a  +  Vo )  * >  0  • 


(6.122) 


so  that  as  t  -*■  oo  , 


q(t)  -v 


nr- n+1 

lL_t-r 


f-oo  (‘n+1 1 


constant 


for  0  <  n  <  1  , 


for  n  =  1  , 


for  n  >  1  . 


(6.123) 


Case  IV:  w"(0)  <  E  ,  crn(t)  -*■  oo 


Since  In  this  case  aQ(t)  eventually  exceeds  the  results  of  the 
analysis  for  Case  III  allow  us  to  conclude  that  the  phase  boundary  will 
become  the  leading  disturbance  at  some  time,  again  denoted  by  t^.  For 
t  >  t^,  the  governing  equations  are  once  more  given  exactly  by  (6.34). 
An  argument  like  that  given  In  Case  I  yields 


■4. 


(6.124) 


This  result  could  have  been  anticipated  by  formally  letting  o  -*■  a>  in 
(6.107).  We  also  find 


f-oo  ■  *<0  • 


(6.125) 


As  In  Case  I,  It  follows  that 


i 
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11m  q(t)  *  +co , 
t-*-® 


(6.126] 


j  a0(>/Hd8'*  -  I^Dt  >  0  *  (6-127> 

Let  us  again  consider  the  asymptotic  power-law  loading  of  Example  1, 
whence 

o  (t)  *  tn  +  o  ( tn )  ,  k.j  >  0  ,  0  <  n  <  a>. 

Then  (6.127)  and  (6.77)  yield 

q(t)  *  (-  ^  D)"71"  ^  t"*1  .  (6.128) 

The  simple  form  of  (6.34)  allows  the  computation  of  the  next  cor¬ 
rections  for  s(t)  and  f(z),  provided  we  know  the  next  term  In  an 
asymptotic  series  for  oQ(t).  Consider  the  following  extension  of  (6.77) 

aQ(t)  *  k.jtn  +  k2tm  +  o(tm)  ,  ^>0,  n  >  0  ,  m<n.  (6.129) 

An  asymptotic  analysis  of  (6.34)  then  reveals  that 

r  1  r  1 

as  t-oo,  (6.130) 

-n 

f(2h^Tl(,?  as  z  -cd  ,  (6.131) 

where 

1  See  the  discussion  leading  to  (6.74). 

2  See  (6.76). 
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7.  Suninary  and  Conclusions 

In  this  final  section  we  shall  describe  the  strain,  stress  and 
velocity  fields  In  the  bar,  with  particular  reference  to  their 
asynptotlc  behavior  for  times  both  very  near  tg  and  much  greater  than 
tg  .  Also  we  present  a  summary  of  our  results  for  the  bar  loaded  In  a 
hard  device.  We  close  with  a  discussion  of  other  pertinent  aspects  of 
the  problems  we  have  considered. 

A.  The  Strains  and  Velocities  on  the  Phase  Boundary 
We  have  already  seen  that 

^G(s(t)",t)  <  0,  v(s(t)",t)  <  0,  t  >  tg,  (7. 1)1 

hence  the  strain,  stress  and  velocity  Immediately  ahead  of  the  phase 
boundary  are  monotonlcally  decreasing  with  time.  Inwedlately  behind  the 
phase  boundary  we  have  from  (4.10)  that 

3t  v(s(t)\t)  .  Vf  -  Vf  S >' 

(7.2) 

-y[Z)  +  C0(t  +  S(t)). 

and 

^€(s(t)tt)  -  -  f'(s(t)  .tftMt)  -yjf)  *  f  (-s(t)  t) 

*  (-S(t)  -,/f)  ♦  E0(t  S(t))(l  S(t)).  7'3) 

From  (2.32),  (6.2),  (6.17),  (6.60)  and  (7.2)  It  follows  that 


See  (6.7). 


1 


v(s(t)  ,t)  >  0  ,  t  >  tg. 


(7.4) 


To  obtain  a  corresponding  monotonicity  result  for  e(s(t)+,t)  we  must  go 
beyond  (7.3)  since  there  the  first  term  works  against  us.  Substituting  from 
(6.53)  and  (6.54)  Into  (7.3)  we  find  after  some  algebra  that 


♦  Jill 


r  L  I  -  ,  J 


2^E  [(s(t)  -Jf)2  (i(t) 


3t  e0(T(s<t)»t)) 


i=~  (w'(c„(x(s(t),t)))  - 


(7.5) 


(7.6) 


Ee0(x(s(t) .t))  -  D). 

Hence  (6.1),  (6.7)  (6.8),  (6.57)  and  (7.5)  yield 
~  e(s(t)+,t)  >  0  ,  t  >  tg  . 

Thus  the  strain,  stress  and  velocity  Immediately  behind  the  phase 


boundary  are  monotonlcally  Increasing  with  time. 


B.  Summary  of  Results  for  Small  Time 


From  (5.40)  -(5.42),  the  small-time  asymptotic  results  for  s(t)  and 
f(z)  are  summarized  as  follows: 


s(t> .  -  (  °°(V  V/z  (t  - 1/'*  .o[(t-  t.m  . 

\3p(6z  -  87)y  6  L  8  J 

f  (i)  - 1  gf  *(8,  >  -  b2>  ♦  *  (f>-3/4  J2 


(7.7) 


(7.8) 
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An  appeal  to  (4.10),  (4.2),  (7.7)  and  (7.8)  determines  the  strain,  stress 
and  velocity  Imnedlately  behind  the  phase  boundary: 


e(s(t)+,t)  *  S2  +  T 


3  ^(‘o) 


(t  -  tg)  +  o[(t  -  tg)]. 


(7.9) 


a(s(t)+,t)  3  ag  +f<70(tg)(t  -  tg)+o[(t  -  tg)]. 


(7.10) 


+  1  /^(Bp-Bl  )<j_(t«)  \l/2  ,  jn  lyo 

r(s(t)+,t)  «  +(0,)  p  P  )  (t-tg)1/2+o[(t-tg)1/2]. 


(7.11) 


As  required,  these  three  field  quantities  are  smoothly  Increasing  with 
time.  In  particular  we  note  that  although  the  quantity  B2-B-|  —  which 
goes  to  zero  with  the  length  of  the  Interval  for  the  unstable  material 
phase  -  appears  In  the  denominator  of  sQ,  It  does  not  occur  In  the 
denominators  of  any  of  the  expansion  coefficients  appearing  In 
(7.9M7.11). 

For  comparison,  we  display  the  leading  terms  In  the  series 
expansions  for  the  prescribed  load  and  the  resulting  strain  and  velocity 


at  the  end  of  the  bar: 

oAta)  o 

e(0,t)  «  B2  +  -2^-2-  (t  -  tg)  +  oC(t  -  tg)2], 


(7.12) 

(7.13) 


o(0,t)  *  Og  +  aQ( tg) (t  -  tg)  +  OC(t-tg)2],  (7.13) 

i  3(Bo“Bi )d  ( t„)  1/2  i/o 

v(0,t)  *  <fr(B-| )  +  ^  - -  1  °  P  (t-  tgr2+0[(t-tg)1/2].  (7.14) 


Between  the  end  of  the  bar  and  the  phase  boundary,  where 


s(t)  <  x  <  0  , 


(7.15) 


we  have  that 


(7.16) 


B2  +  T-2^"  (t"t8)  *  oC(t>t8^  -  e(x,t)  -  82  +  ~2TL~  (t‘t8) 


+  o[(t-t6)]. 


°8  +  t  W^"V  +  -  0^x»t^  i  °e  +  ^V^V' 


+  o[(t-tg)]. 


(7.17) 


v(x,t)  «  <j>(8-|)  +  j 


i  /acs^JooCt  }  \V2 


)1/2(t-te)1/2  +  o[(t-te)1/Z3.  (7.18) 


We  note  that  the  proposed  expansion  scheme  (5.11),  (5.12)  and 
(5.43)  would  suggest  that  the  Ho"  error  bounds  appearing  In  (7.7)  -<7,11 ), 
(7.16)-(7.18)  may  be  tightened  to  "O"  of  the  next  larger  half-integer 


poMer. 


We  turn  now  to  a  determination  of  the  various  field  quantities 
Immediately  ahead  of  the  phase  boundary.  Here  from  (5.22),  (5.23),  (5.39) 


and  (5.40): 


'  1;^, *  •««A 


(7.19) 


Furthermore  (2.30),  (5.25),  (5.26),  (5. 29) -(5. 31 ) ,  (5.39)  and  (5.40)  yield 
e(s(t)",t)  *  6-,  -( - 2— - )(t-tfi)  +  oC(t-t-)],  (7.20) 

1  \2(6Z^  W"(^))  6  8 

o(s(t)-,t)  -  o„  -  (■  -Oo-V„ - Yt-tJ2  +  o[(t-tJ2],  (7.21) 

8  \  18(82^  rw  (B1)/  8  8 

\  _  x/n  \  (  '  2  °o(t8)3/2  \  *  \3/2  A 


v(s(t)',t)  -  4(8,)  -  ^-T/2("°.J/rw,„(6i)j  («■>"  + 

+  o[(t-tJ3/2].  (7.22) 


As  required,  the  above  field  quantities  are  smoothly  decreasing  with 
time.  He  also  note  the  presence  of  the  denominators  of  the 

coefficients  of  the  terms  giving  the  leading  order  time  dependence.  The 
final  quantity  which  we  shall  examine  Is  the  acoustic  speed  Immediately 
ahead  of  the  phase  boundary.  From  (4.20h),  (7.7)  and  (7,19)  It  follows 


that 


c(e(s(t)',t)) 


W 

3p(82"®i ^ 


[1/2 


(t-t6)1/Z  +  o[(t-t8),/2].  (7.23) 


We  note  that  (7.7)  and  (7.23)  confirm  (6.6)  fbr  small-time.  The  determination 
of  the  exact  order  of  the  remainder  terms  appearing  In  (7.20)-(7.23)  Is 
complicated  by  the  multiple  nesting  of  functional  arguments  occurring 
In  the  various  formulae,  as  well  as  by  the  need  to  determine  the  next 
order  correction  to  (7.19). 

C.  Suninary  of  Results  for  Large  Time 

He  have  treated  separately  the  four  cases  (6.35)-(6.38  )  and  have 
found  that 


where 


and 


s(t)  *  -  at  +  q(t) 

as  t  +  •  , 

(7.24) 

f(z)  *  f^  +  f,(z) 

as  z  +  -0®, 

(7.25) 

o  <  a  < J- 
P 

(7.26) 

q(t)  *  o ( t ) ,  q(t)  =  o(l) 

as  t  +  •  , 

(7.27) 

f-,(z)»  o(T) 

as  z  +  — 

( 7 1 28 ) 

We  first  considered  the  cases  given  by  (6.35)  and  (6.36),  for  which 


For  these  cases  the  phase  boundary  cannot  become  the  leading  disturbance. 
When  the  applied  load  approaches  a  finite  limit  the  asymptotic  phase 
boundary  speed  a  Is  given  by  the  unique  root  of  (6.89)  occurring  within 
the  Interval  o«x<yj .  It  follows  that  a  Is  an  Increasing  function  of 
am  and  approaches  as  o^  tends  to  Infinity.  The  constant  f^  may 
then  be  found  from  either  (6.87)  or  (6.88).  When  the  applied  load  Is 
unbounded,  a  *  -J—  and  11m  q(t)  *  ®.  In  particular  the  dominant 
behavior  of  q(t)  Is  given  Implicitly  by  (6.76),  while  f^  Is  given  by 
(6.73). 

For  the  cases  given  by  (6.37)  and  (6.38),  we  have 

WH(0)  <  E.  (7.30) 


For  these  cases.  If  the  applied  load  remains  below  the  value  o^,  given 
by  (6.110),  the  phase  boundary  does  not  become  the  leading  disturbance. 
Here  a  Is  the  unique  root  of  (6.89)  occurring  within  the  Interval 
(0,c(0)],  while  f^  Is  again  found  from  either  (6.87)  or  (6.88).  If  the 
applied  load  exceeds  o^  the  phase  boundary  eventually  becomes  the  leading 
disturbance.  In  particular  If  the  applied  load  approaches  a  limit 
a  >  c, ,  a  Is  given  by  (6.107)  and  f  by  (6.112).  If  the  applied  load  is 
unbounded,  a  ■/-■  and  f^  *  0/4E.  In  this  latter  circumstance  q(t) 

Is  also  unbounded,  the  dominant  behavior  being  given  implicitly  by 
(6.127). 

We  have  already  noted  that  ahead  of  the  phase  boundary, 
e(s(t)",t),  o(s(t)",t)  and  v(s(t)",t)  are  monotonlcally  decreasing  with 
time.  Of  course  If  the  phase  boundary  becomes  the  leading  disturbance, 
say  again  at  ta  t^,  these  quantities  are  all  zero  for  t  >  t^.  Otherwise 
the  phase  boundary  Is  confined  to  the  disturbed  region  and  is  asymptoti- 
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Conslder  the  fixed  particle  P  Initially  at  x.  Since  s(tQ)  *  0,  It 

D 

follows  from  (5.9)  and  (6.1)  that  for  any  given  x  <  0  the  equation 
x  *  s(t)  has  a  unique  root  t  ■  tx  with  tx  >  tg.  Consequently  the 
particle  P  finds  Itself  In  the  high-strain  phase  for  times  greater  than 


tx.  For  fixed  x  we  have 


Urn  f(x  ->/|  t  )  -  f^  ,  11m  f(-x  -V|  t 


(7.35) 


t-*»  ~  t-»- 

whlch,  with  (4.6)  and  (4.10),  yield  the  following  large-time  asymptotic 
results  for  the  stress,  strain  and  velocity  at  the  particle  P  initially 
at  x: 


It  Is  Interesting  to  compare  (7.36)-(7.38)  to  the  stress,  strain 
and  velocity  fields  occurring  In  a  companion  problem  Involving  a 
linearly  elastic  bar,  for  which  the  stress  Is  related  to  the  strain  via 
1 - 

See  (6.11),  (6.12)  and  (6.28). 


cally  parallel  to  the  characteristic  ray  C_  In  Rj  which 

x  ■  0  at  time  t  *  r  ,  where 

Issues  from 

.  -  V  •.*« 

T.  *  fl(ci)  . 

Then  (6.30)-(6.32)  yield  the  asymptotic  results: 

(7.31  )7 

c(s(t)‘,t)  <\i  c"7(a)  as  t  ®, 

a(s(t)~,t)  ^  W (c-1  (a))  as  t  ■*■  «, 

(7.32) 

-  ■  -Vj 

(7.33) 

v(s(t)",t)  <fr(c"7(a))  as  t  ®. 

(7.34) 

-*.*•  „** ' 

a(x,t)  -v/  aQ(t  +Vf  x) 

as  t  ®, 

(7.36) 

e(x,t)  -v  £Q(t  +Vf  x) 

*  [a0(t+VFx>  "  D^/E 

as  t  ®, 

(7.37) 

v(x,t)  [aQ(t+Vf  x) 

as  t  ®. 

(7.38) 

UiA 

..  *  ■ 
.  • .  * . 


I 


o  3  Ee  +  Dg.  e  >.  0. 


(7.39) 


In  this  companion  problem  we  shall  denote  the  strain,  stress  and  velocity 
fields  by  e(x,t),  o(x,t)  and  v(x,t).  The  response  of  such  a  linearly 


elastic  bar  subject  to  the  prescribed  end-loading 


3(0,t)  3  cr0(t). 


In  the  disturbed  region  t  is  give 

3(x,t)  3  aQ(t  +  |  x)  , 

e(x,t)  3  [aQ(t  +^|  x)  -  D2]  /  E  , 

v(x,t)  3  CaQ(t  +^f  x)  -  D2]  pE  . 


ven  by 


(7.40) 


(7.41) 


(7.42) 


(7.43) 


Two  cases  of  Interest  In  comparing  (7.36>*l7.38)  to  (7.41  )-(7 .43)  are 
D2  3  0  -  In  which  case  the  linearly  elastic  bar  when  stress-free  Is 
unstrained  -  and  D2  3  D,  in  which  instance  ( 7. 39)1  s  the  linear  extrapolation 
back  through  e»0  of  the  final  linear  portion  of  the  stress-strain  curve 
occurring  In  the  original  problem.  Recall  that  D  can  be  either  positive 


or  negative  ,  According  to  (6.63), 

so  that  tg)  may  be  either  positive  or  negative  depending 

on  the  value  of  02  and  on  the  detailed  behavior  of  W(e)  within  the 
interval  (0,B-|).  Similarly  f^  may  be  either  positive  or  negative2 

As  we  have  seen  D<0  when  W"(0)<E;  see  (6.97). 


(7.44) 


2  For  example  when  W" ( 0)<E,  1 1  follows  that  f(-^j^  tg)<0.  However  If 
W"(e)>9E  for  0<e<-jr  6-| ,  and  e2  <  \  B-j ,  it  follows  that  f(-  ts)>°. 

In  this  case  (7.44)  provides  f^O.  Conversely, (6.112)  furnishes  an 
example  with  f  <0. 
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Con  sequent!  y,  apart  from  the  monotonlclty  of  f,  there  Is  little 
Me  can  say  concerning  Its  range  without  stipulating  further  conditions 
on  the  stress-strain  relation. 

One  case  of  Interest  Is  that  of  a  material  with  a  small 
nonlinearity.  Suppose 

|W"(e)  -  E|  £  Sg  f°r  0  <  e<Bi-61#  and  that  82  a  8-j  +  5-j»  (7.45) 

where 

B-j  »  6^  >  0  and  E  »  62  >  (7.46)1 

Here  the  material  Is  "almost  linear"  except  In  a  S-j -neighborhood  about 
e*  6-j.  For  such  a  material 

11m  f(-J|t8)*  11m  i  //6l^H5d8-B\  =0,  (7.47) 

6r62-0  6  5r62-K)Z  \0  E  Z/ 

/ 

and 

11m  D  *  11m  (Wl(g, )  -  E6o)  =  0  ,  (7.48) 

5-j  ,62-*0  6^»52-K) 

We  now  show  that  In  all  of  the  cases  set  down  in  (6.35)-(6.38) , 

11m  f^  *  0  .  (7.49) 

5-j  ,<$2-*0 

When  the  phase  boundary  becomes  the  leading  disturbance, (7. 49)  follows 
from  (7.48)  and  either  (6.112)  or  (6.125).  When  the  phase  boundary  is 
confined  to  the  disturbed  region  ,  we  examine  either  (6.73)  or  (6.88). 
First  note  that  for  e  <  8, 

1  If  It  were  not  for  (6.18),  we  would  take  <S2=0. 
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lim  W/(e)  *  Ee  , 

6-|  ,52-*^ 


(7.50) 


11m  4»(e)  *  lim  /  -J  ^  ds  =  /  ds  *  S^-  e.  (7.51 ) 

.  r*  m  •.  n  />*  ”  k  n  "  0  "  Q 


6-j  »52‘^0 


51  »<s2'*0  a 


When  (6.88)  applies*  we  have 
11m  f^  *  —i 

&1 


71 *,1^0  ft  •W!“c‘1<“)+7<“+a/|)*(c'1<“)) 

-  fe^'wi] 

+  ^6  1J\0(^ f<KC'1(a))  "  P  y/  (c’1(a)))]=  °  •  (7-! 


Similarly  (7.48),  (7.50)  and  (7.51)  yield  (7.49)  when  (6.73)  applies. 
This  completes  the  verification  of  (7.49)  In  all  cases.  Finally  from 
(7.44), (7,47)  and  (7.49)  we  have 


lim  f(z)  =  0 
<$1 ,62'^ 


(7.53) 


It  is  Immediate  from  (7.45)  that'  the  strain,  stress  and 
velocity  ahead  of  the  phase  boundary  will,  as  6^  and  62  approach  zero, 
approach  the  corresponding  values  occurring  In  a  material  with  the 
stress- strain  relation 

05  Ee  e  >  0  .  (7.54) 


What  (7.53)  guarantees  Is  that  the  same  Is  true  behind  the  phase 
boundary^  .  Thus  It  follows  that  the  jumps  In  strain,  stress  and 
velocity  across  the  phase  boundary  win  tend  to  zero  with  6-j  and 
In  other  words  the  solution  to  the  nonlinear  problem  approaches  the 
“shock-free"  solution  to  the  linear  problem  as  the  nonlinearity  becomes 
small2. 

D.  Sunmary  of  Results  for  Hard  Loading 

So  far  we  have  only  considered  "soft"  loading  conditions  for  the 
materials  Introduced  In  the  fourth  section.  For  hard-loading  ,  we  have 
formulated  In  (3,24)  the  free  boundary  problem  governing  the  phase 
boundary  for  materials  unrestricted  by  assumptions  as  to  the  curvature 
of  the  stress-strain  relation  In  the  high-strain  material  phase.  As  one. 
would  expect,  tfwe  con  fine  our  attention  to  the  subclass  of  materials  for 
which  (4.2)  also  holds,  we  are  able  to  Integrate  (3.24)  partially  whereupon 
we  arrive  at  a  pair  of  differential -delay  equations  formally  resembling 
(4.22)  and  (4,23),  We  may  then  recast  the  problem  as  an  Integral 
equation  formally  resembling  (4.32). 

Proceeding  as  In  the  analysis  of  the  soft,  device,  we  arrive  at 
qualitatively  similar  results.  The  major  difference  between  hard  and 
soft  loading  occurs  In  the  small-time  asymptotic  results.  We  restrict 
attention  to  the  generic  case:  W^S-i)  <  0  and  v  (tfl)  >  0.  Here  t  *  t„ 

I  Op  p 

Is  the  unique  root  of  (3,23).  Upon  completion  of  an  analysis  parallel  to 
that  of  the  fifth  section,  one  Is  led  to  the  following  asymptotic  results: 

T - 

Formally  set  f»0  In  (4.10), 

2 

This  Is  reassuring,  since  naively  letting  in  (7.7)  and  (7.20)- 

(7,22)would  have  led  us  astray, 


S(t) 


(7.55 


»  -  t8r  +oi(t  -  tgn  , 

2(02-0!)  B  6 

.  pV2(tfl)  9  o 

e(s(t)+,t)  *  82  + - — (t  -  tft)Z  +  o[(t  tft)Z]  \  (7.56] 

2  2(02-0!)  8  6 

p.2(t 

a(s(t)+*t)  -  a.  +  -9J&L  (t  -  t  )2  +  o[(t  -  t  )2j  ,  (7.57] 

p  62-6-| 

v(s(t)+,t)  «  v0(tg)  +  vQ(tg)(t  -  tg)  +  0[(t  -  tg)2]  .  (7.58] 


Also,  In  analogy  with  (7.16)  arid  (7.18),  here  we  have  for 


that 


x  <  s(t)  <  0 


(7.59] 


Pv?(tft)  9  9  3pv2(tJ 

6  (t*tj2  +  o[(t-tfl)2]<  e(x,t)  <  B2  + — g— 1 — 

2  £(62-6^  B  8  2 


2E(62"6] ) 


(t-t-)‘ 


+  oC(t-  tg)2] 


(7.60] 


and 

v(x,t,)  *  vQ(t)  +  o[(t-tg)2]  .  (7.6V 

Moreover,  In  contrast  to  the  soft  loading  case,  Indications  are  that  here  the 
phase  boundary  location  s(t)  —  as  well  as  the  strain,  stress  and 
velocity  fields  Immediately  on  either  side  of  the  phase  boundary— can 
apparently  be  expanded  In  Integer  powers  of  t-tg. 

Turning  to  the  large-time  analysis  we  again  represent  s(t)  In  the 
form  (7,24)  subject  to  (7.26)  and  (7.27).  As  with  soft-loading,  the 
phase  boundary  cannot  become  the  leading  disturbance  for  materials  with 
W“(0)  >  E,  In  this  case,  If  the  end  velocity  Is  unbounded,  a  ‘yfjf  and 
q(t)  Is  unbounded,  while  If  v  (t)  approaches  a  limit  vm  we  find  a  as  the 
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unique  root  within  the  Interval 


(o4)°f 


the  equation 


(a2-  £)rva,-*(c'1(a))]  +  2-[W'(c'1(ct))  -  Ec’1  (a)  -  D]  .  0.  (7.62) 


For  materials  with  Ww(0)  <  E  the  phase  boundary  becomes  the  leading 
disturbance  If  and  only  If  the  Increasing  velocity  at  the  end  of  the 
bar  eventually  exceeds  the  value  v»v^,  where 


(7.63) 


If  vQ(t)  Is  unbounded,  a  and  q(t)  Is  unbounded?  otherwise  a  Is 
found  from  (7,62)  when  v^  <  vL,  while  It  Is  given  by 


when  >  vL  ,  In  all  cases  It  Is  possible  to  obtain  higher  order 
corrections  through  methods  analogous  to  those  employed  In  Section  6. 


E.  Further  Remarks 

The  Integral  equation  formulation  for  this  problem  remains  relatively 
unexplored.  It  can  be  shown  that  for  a  certain  type  of  plecewise-1  Inear 
stress-strain  relation,  the  Integrals  can  be  evaluated  yielding  a  scalar 
equation  relating  the  quantities  t,  tfl  and  tb.  Such  an  equation  could 
perhaps  permit  the  determination  of  s(t).  It  should  be  noted  that  If 
s(t)  Is  known,  It  Is  a  simple  matter  to  determine  u(x,t)  In  by  &  . 
procedure  we  now  briefly  outline.  Since  (4.15)  ensures  that  the  curve 
x  *  s(t)  Is  never  parallel  to  the  characteristic  rays  In  t^,  the  phase 
boundary  Is  an  acceptable  curve  on  which  to  specify  Cauchy  data  for  the 
classical  wave  equation  (4.11)  which  governs  u(x,t)  In  Rg.  The  proper 


Cauchy  data  to  specify  on  the  phase  boundary  may  be  found  from 
e(s(t)+,t)  and  v(s(t)+,t), which  In  turn  are  given  In  terms  of  the  known 
functions  e(s(t)",t),  v(s(t)",t)  and  s(t)  by  means  of  the  jump  conditions 
(2.12).  Then,  as  Is  well  known,  an  application  of  Green's  theorem  In  the 
plane  will  yield  u(x,t)  In  It  Is  essential  to  locate  the  curve 

x  *  s(t)  properly  In  order  to  guarantee  that  this  solution  to  the  classical 
wave  equation  will  satisfy  the  condition  prescribed  on  x  *  0. 

If  this  procedure  Is  employed  with  s(t)  as  yet  unknown.  It  gives 
a  functional  representation  for  u(x,t)  behind  the  phase  boundary 
In  terms  of  s(t)  and  s(t),  For  the  problem  of  hard  loading,  equating 
the  prescribed  value  of  u(0,t)  to  that  provided  by  this  functional 
representation  furnishes  an  elegant  alternative  derivation  of  the  Integral 
equation.  However,  since  the  soft  boundary  condition  amounts  to  prescribing 
the  normal  derivative  of  u(x,t)  on  x  ■  0,  this  method  will  not  supply 
the  Integral  equation  for  the  case  of  soft  loading. 

One  notable  aspect  of  our  analysis  Is  the  Inequality  (6.8)  which 
assures  that  the  phase  boundary  Is  subsonic  with  respect  to  the  material 
In  the  high-strain  phase  behind  It,  but  supersonic  with  respect  to  the 
material  In  the  low-strain  phase  ahead  of  It.  According  to  James  [19] 
certain  uniaxial  tensile  experiments  on  a  bar  In  a  dead-loading  (soft) 
device  Indicate  that  phase  boundaries  appear  which  separate  approximately 
homogeneously  deformed  phases.  These  experiments  also  indicate  that  the  phase 
boundary  travels  more  slowly  than  the  sound  speed  of  either  phase.  That  the  phase 
boundary  we  describe  Is  supersonic  with  respect  to  the  material  ahead  of 
It  Is  equivalent  to  the  geometrical  requirement  that  the  leading  side 
of  the  phase  boundary  Intersect  characteristic  rays  which  arise  either 
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at  the  end  of  the  bar  x«0  for  t<tg,  or  on  the  line  of  Initial  data  x<0,  t*0. 
From  the  discussion  following  (3.5)  we  recall  that  the  vanishing  of  the 
sound  speed  at  x  *  0,  t  *  tg  ensures  that  these  characteristic  rays  form 
an  unambiguous  cover  of  (0, t) 1 1  >  tg}.  Consequently  any  phase 
boundary  emerging  at  time  tg  must  at  least  Initially  be  supersonic  with 
respect  to  the  material  ahead  of  1t^  .  Moreover,  when  (6.1)  holds  the 
curvature  of  the  phase  boundary  In  the  x-t  plane  assures  that  the  leading 
side  of  the  phase  boundary  continues  to  Intersect  new  members  of  this 
covering  set  of  characteristic  rays,  thus  ensuring  the  supersonic 
condition  (6,6), 

For  a  phase  boundary  to  advance  subsonlcally  Into  a  material  In  the 
low-strain  phase,  characteristic  curves  must  arise  on  the  leading  side 
of  the  phase  boundary.  Such  a  violation  of  the  entropy  condition  given 
by  Lax  would  seem  possible  If  the  phase  boundary  were  to  emerge  on  x  *  0 
before  the  sound  speed  of  the  material  vanishes,  l.e.  before  the  loss 
of  hyperbollctty  at  time  tg.  Suppose  the  phase  boundary  emerges  on  x  =  0 
at  time  t0  <  tg.  Then  a  sufficient  condition  for  the  phase  boundary 
to  be  Initially  subsonic  with  respect  to  both  material  phases  Is  that 
It  emerge  with  zero  velocity.  In  this  circumstance  any  point  In  the 
region  of  the  x-t  plane  between  the  characteristic  ray  originating  at 
x  ■  0,  t  *  t  and  the  phase  boundary  would  lie  on  a  negatively-sloped 
characteristic  ray  which  originates  on  the  phase  boundary.  Consequently 
the  solution  representation  (2.30)  would  apply  only  up  to  the  character¬ 
istic  ray  originating  at  x  *  0,  t  *  te.  The  obvious  candidate  for  a 
solution  In  the  region  -c(e0(te))(t-te)  <  x  <  s(t)  would  be  the  extension 
1 

As  In  confirmed  by  (7.7)  and  (7.23). 
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of  the  strain  and  velocity  fields  as  constants  from  x*  -c(eQ)(te) ) (t-tg) . 

The  time  at  which  the  phase  boundary  emerges  Is  determined  by  the 
particular  Inverse  to  the  stress  response  function  a*  Wy (e)  which 
operates  at  the  end  of  the  bar,  In  our  study  we  have  assumed  an  Inverse 
which  results  In  aehangeof  phase  only  when  the  current  phase  can  no 
longer  support  the  applied  load.  A  consequence  of  this  assumption  Is  that 
the  sound  speed  of  the  original  phase  vanishes  where  the  phase  boundary 
emerges.  The  vanishing  of  the  sound  speed  Is  a  phenomenon  not  without 
precedent  In  continuum  mechanics.  Consider  the  Isentroplc  one-dimensional 
non-steady  flow  engendered  by  withdrawing  a  piston  from  a  tube  containing 
an  Ideal  polytropic  gas  Initially  at  rest.  If  the  final  piston  velocity 
exceeds  the  escape  speed,  the  sound  speed  Is  zero  at  the  leading  edge  of 
the  gas  as  it  moves  into  the  newly  created  volume.  Although  this  piston 
problem  has  many  formal  similarities  to  our  problem  of  the  elastic  bar, 
there  is  no  appropriate  analogue  of  a  second  phase.  Here  the  region 
between  the  piston  and  the  leading  edge  of  the  gas  Is  simply  a  vacuum  or 
zone  of  cavitation.  A  relevant  solid  mechanical  analogue  of  the  piston 
problem  would  seem  to  be  the  monotonic  loading  of  an  elastic  bar  whose 
stress  response  curve  Is  Initially  Increasing  with  strain  and  then 
decreasing  without  a  final  Increasing  portion,  for  then,  just  as  speeds 
above  the  escape  speed  cannot  be  achieved  by  the  gas,  equilibrium 
stresses  above  a  certain  level  cannot  be  sustained  by  the  elastic 
solid1  , 

1  For  a  discussion  of  such  a  stress  response  In  (Incompressible)  elastic 
solids  and  the  corresponding  analogy  with  compressible  fluids,  reference 
may  be  made  to  Section  4.  of  [23]. 
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As  discussed  In  the  third  section.  It  Is  possible  to  Introduce  a 
second  phase  emerging  at  the  end  of  the  bar  at  any  time  In  the  Interval 
[t  ♦to]1  .  Besides  our  choice  of  t  *  tfl,  another  plausible  candldsce  Is 

<*P  P 

the  time  t^  when  the  applied  load  reaches  the  Maxwell  strts*.  Then 
o0(tn|)  *  oy  and  we  expect  to  observe  o  *  ay,  e  *  y-j,  and  v  *  4>(y1)  In 
the  region  of  the  x-t  plane  between  the  "Maxwell  characteristic"  and 
the  phase  boundary.  Although  we  have  yet  to  Investigate  such  a  problem, 
we  point  out  that  since  the  fields  In  front  of  the  phase  boundary  would, 
at  least  Initially,  be  homogeneous,  such  a  study  has  the  analytical 
advantage  of  eliminating  the  need  to  Invert  (4.20h), 

Besides  the  Maxwell  characteristic,  there  Is  another  canonical 
characteristic.  We  observe  that,  at  the  end  of  the  bar,  the  particles 
are  Initially  subsonic  since  v(0,0)  *  0  <  *  c ( 0, 0) ;  at  time  tg, 

however* they  are  supersonic  since  v(0,tg)  *  4(8^  >  0  *  c(0,tg).  Hence 
there  Is  a  “transonic  characteristic"  on  which  the  particles  are  moving 
to  the  right  at  the  acoustic  speed.  We  can  locate  this  characteristic 
by  solving  for  the  associated  "transonic  strain".  The  latter  Is  the  unique 
root  within  the  Interval  (0,8^)  of  the  equation 


c(e)  -  <p(e)  *  ^W^(c)  -  /  ^W"(e)  de^*  0. 


(7.65) 


In  contrast  to  the  familiar  situation  of  two-dimensional  Isentroplc  steady 
compressible  fluid  flow,  the  governing  equations  for  the  dynamics  of  an 
elastic  bar  do  not  change  their  type  as  the  continuum  passes  from 
subsonic  to  supersonic.  It  should  be  noted  from  (7.65)  that  the  transonic 


See  the  discussion  surrounding  eqn.  (3.8). 
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strain  Is  completely  determined  by  the  first  ascending  branch  of  the 
stress-strain  curve.  On  the  other  hand,  the  Maxwell  strain  for  the 
material  In  the  low-strain  phase  depends  on  the  global  stress-strain 
relation  via  the  equal  area  rule.  Consequently  the  transonic  strain  may 
be  either  greater  than,  equal  to,  or  less  than  Yj.  In  other  words, 
whether  the  particles  become  supersonic  before  or  after  they  reach  the 
Maxwell  stress  depends  on  the  particular  stress-strain  curve  under 
consideration. 

Another  Interesting  result  of  our  analysis  Is  that  the  velocity  of 
the  phase  boundary  approaches  the  sound  speed  of  the  high-strain  phase 
^  as  the  applied  load  tends  to  Infinity.  This  bears  a  resemblance 
to  results  In  other  free-boundary  problems  In  mechanics.  Burrldge 
and  Keller  [24]  treat  a  semi-infinite  one-dimensional  linear-elastic 
adhesive-tape  which  Is  Initially  completely  attached  to  a  plane  surface 
and  subject  to  an  In-plane  tensile  force  at  Its  end.  Peeling  Is  Initiated 
by  suddenly  applying  a  sufficiently  large  normal  force  to  the  end  of 
the  tape.  The  free-boundary  Is  the  point  separating  the  portion  of  the 
tape  still  adhering  to  the  surface  from  that  part  which  has  already 
peeled.  If  the  pair  of  forces  at  the  end  of  the  tape  is  maintained,  the 
free-boundary  will  move  at  a  constant  velocity  which  Is  a  monotonically 
Increasing  function  of  the  magnitude  of  the  applied  normal  force. 

Moreover  the  value  of  this  function  approaches  the  fixed  sound  speed  of 
the  elastic  adhesive-tape  as  the  magnitude  of  the  applied  normal  force 
tends  to  Infinity.  The  analogy,  however,  cannot  be  carried  too  far,  since 
the  tape  problem  incorporates  a  condition  which  characterizes  the 
breaking  of  the  adhesive  bond  between  the  tape  and  the  surface,  while  the 
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problem  of  the  elastic  two-phase  bar  Incorporates  no  analogous  failure 
condition  external  to  elasticity  theory.  The  adhesive  failure  condition, 
which  applies  only  at  the  free-boundary.  Is  dissipative  In  the  sense  that 
It  lowers  the  free-boundary  speed^.  It  would  seem  that  a  corresponding 
effect  on  the  phase  boundary  speed  might  occur  In  the  problem  treated 
here  If  a  surface  energy  were  associated  with  the  phase  boundary.  As 
James  [5]  has  observed,  the  effects  of  such  surface  energies  have  so 
far  been  Ignored. 

Finally  we  note  that  In  this  study  we  have  totally  neglected  unloading, 
a  process  which  we  would  expect  to  promote  conventional  shocks;  see  [16]. 

In  addition  we  have  dealt  with  a  very  specialized  class  of  materials  that 
not  only  made  possible  the  Integration  of  some  of  the  governing  equations, 
but  also  precluded  conventional  shocks  behind  the  phase  boundary.  For 
more  general  materials,  even  determining  the  smoothness  of  the  fields 
behind  the  phase  boundary  would  perhaps  demand  mathematical  techniques 
entirely  different  from  those  employed  In  this  study;  see,  for  example,  [25]. 
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Figure  2.  Geometry  of  the  characteristic  rays  C_  for  an  elastic  singl 
phase  bar. 


Figure  3a.  Stress-strain  curve  with  an  Isolated  discontinuity  in  W"(e) 
at  e  =  e 


C_(td+)  :  x  =  -c(«d+)(t- td) 
C_(td“)  :  x  =  -c(€d")  (t-  td) 
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Figure  3b.  Region  between  the  characteristic  rays  associated  with  strai 
ed  which  involves  constant  strain  and  velocity. 
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Figure  5.  The  smooth  monotone  loading  past  o*og  (*o0(tg))  of  a  bar 
whose  material  is  non-hardening  in  the  first  phase  is  as¬ 
sociated  with  a  pattern  of  characteristic  rays  C_  which  pro 
vide  an  unambiguous  covering  of  the  region  fc\l(0,t)|t*tQj. 


Figure  6.  The  particular  Inverse  to  the  stress  response  curve  of  Figure  1 
which  results  in  the  material  assuming  the  high-strain  phase 
only  when  the  low-strain  phase  can  no  longer  support  the  ap¬ 
plied  load. 
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whose  stress-strain  relation  Is  not  monotonic.  Sufficiently  large  applied 
loads  then  require  the  strain  to  jump  from  one  ascending  branch  of  the  stress- 
strain  curve  to  another  such  branch.  For  a  special  class  of  these  materials, 
a  nonlinear  Initial-boundary  value  problem  In  one-dimensional  elasticity  Is 
considered  for  a  semi-infinite  bar  whose  end  Is  subjected  to  either  a  monotonl- 
cally  Increasing  prescribed  traction  or  a  monotonlcally  Increasing  prescribed 
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displacement.  If  the  stress  at  the  end  of  the  bar  exceeds  the  value  of  the 
stress  at  any  local  maximum  of  the  stress-strain  curve  a  strain  discontinuity 
or  "phase  boundary"  emerges  at  the  end  of  the  bar  and  subsequently  propagates 
Into  the  Interior.  For  classically  smooth  solutions  away  from  the  phase 
boundary,  the  problem  Is  reducible  to  a  pair  of  differential -del ay  equations 
for  two  unknown  functions  of  a  single  variable.  The  first  of  these  two  functlc 
gives  the  location  of  the  phase  boundary,  while  the  second  characterizes  the 
dynamical  fields  In  the  high-strain  phase  of  the  material.  In  these  equations 
the  former  function  occurs  In  the  argument  of  the  latter,  so  that  the  delays 
In  the  functional  equations  are  unknown.  A  short-time  analysis  of  this  system 
provides  an  asymptotic  description  of  the  emergence  and  Initial  propagation 
of  the  phase  boundary.  For  large-times,  a  different  analysis  Indicates  that 
the  phase  boundary  velocity  approaches  a  constant  which  depends  on  material 
properties  and  on  the  ultimate  level  reached  by  the  applied  load  as  well. 
Higher  order  corrections  depend  on  the  detailed  way  In  which  the  load  Is  applie 
Estimates  for  the  various  dynamical  field  quantities  are  glv^n  and  a  priori 
conditions  are  deduced  which  determine  whether  the  phase  boundary  eventually 
becomes  the  leading  disturbance. 
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Figure  8.  Geometry  of  the  characteristic  rays  C_  for  the  material  whose 
stress  response  curve  Is  depicted  In  Figure  7. 
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whose  stress-strain  relation  Is  not  monotonic.  Sufficiently  large  applied 
loads  then  require  the  strain  to  jump  from  one  ascending  branch  of  the  stress- 
strain  curve  to  another  such  branch.  For  a  special  class  of  these  materials, 
a  nonlinear  Initial -boundary  value  problem  In  one-dimensional  elasticity  is 
considered  for  a  semi-infinite  bar  whose  end  Is  subjected  to  either  a  monotoni- 
cally  Increasing  prescribed  traction  or  a  monotonically  Increasing  prescribed 
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